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FOREWORD 


Air Pollution knows no boundaries and as such 
atmospheric pollutants are readily transportable and dis- 
persed by the wayward action of winds. No one group has 
the sole right to contaminate our greatest resource, Air, 
to the detriment of his surrounding neighbors and because 


of this there is an obligation for control at the source, 


Sergius L, Dobko, P. Eng., 
Head, Air Pollution Control Section, 
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AIR POLLUTION IN GENERAL 


INTRODUCTION 


The air that we breathe is one of the most important resources 
we have. Man is directly dependent upon it for breathing and indirectly 
depends on it for his economy and well-being. 


Atmospheric pollution is not a new phenomenon and it undoubtedly 
originated when man first lit a fire. A bush or forest fire must have 
occurred in an area of restricted ventilation resulting in contamination of 
the air by smoke. Air pollution was recognized as a problem during Roman 
times and by the 18th century it was prevalent in most major cities of 
Europe and North America. 


When we think of the air pollution problem, we normally associate 
its source with some activity of man. Man, however, is not solely respon- 
sible for air pollution. For example, the haze normally seen in the country- 
side is associated with plant life. Wind-borne dusts, salt spray and pollens 
are also natural contaminants of the atmosphere. Practically all air is 
contaminated to some extent and air pollution basically is caused by intro- 
duction of foreign substances into the atmosphere. If the concentration 
of these contaminants interferes with the well-being of people and their 
economies, a pollution problem may arise. One of the difficulties in com- 
bating air pollution is the variety of its effects on people. A housewife 
may have her furniture and clothing soiled, a farmer is interested in the 
effect it may have on his crops, a tourist may be inconvenienced by a re- 
duction in visibility or an obnoxious odor. 


Air pollution may be caused by many substances coming from a vari- 
ety of sources. It can come from automobile exhaust pipes, carburetor, 
crankcases, and home furnaces. Waste disposal operations such as incinera- 
tors, open burning of dumps, the burning of leaves and wastes of construc- 
tion and demolition are significant sources of pollution. Large volumes’ 
of pollutants are also emitted from oil refineries, power plants, steel 
mills, pulp and paper factories and almost every other kind of industry. 


Normally, natural atmospheric conditions absorb these pollutants. 
Thorough mixing brings about a dilution of the contaminants; rain, snow, 
and dew helps remove the dust, smoke, and gaseous pollutants. Absorption 
of gases by vegetation, soil, and water are other natural air cleaning 
processes. If this happens as fast as pollutants are released, then the 
consequences are negligible. However, peculiar weather conditions may 
interfere with this normal phenomena and a pollution problem may arise 
due to stagnant air or limited ventilation. 


NEED OF CONTROL 





The population, commerce, and industry in our urban areas is in- 
creasing at a very fast rate and can only be expected to continue to in- 
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crease, Statistics are sometimes felt to be remote and uninteresting, 
but the statistics for the Province of Alberta and the larger urban areas 
of Calgary and Edmonton are very significant to the subject of air pollu- 
tion. 


This rapid expansion in urban areas not only in Alberta but 
throughout the world means that there are a greater number of sources of 
air contaminants which may pollute the atmosphere. This pollution will 
affect us in two ways; firstly, there will be the direct effect on pro- 
perties immediately adjacent, and with the continued development of the 
area, there will be fewer and fewer cases where the adjacent property is 
not developed. Secondly, there will be the cumulative effect of all the 
"sources! within the metropolitan area and as the area becomes larger 
this is a very real consideration, 


During the past few years almost every large urban area in the 
world has experienced a drastic increase in the pollution of its atmosphere. 
As the discomforts of air pollution become more noticeable, public clamor 
for rigid regulation of air contaminating processes increases steadily. 


METHODS OF CONTROL, 





Since it is difficult in generalizing all plants' sources and 
controls under a single report, this section is indicative of the means 
that can be taken of controlling air pollution sources, 


I Fundamental Approaches to Control: 


Regardless of the air pollution problem to be attacked . . that 
is, whether it be a community-wide problem or merely a single-source pro- 
blem . . there are two fundamental approaches to control, the first being 
control of the pollutant at the source so that excessive amounts are not 
emitted to the atmosphere in the first place and the second is natural 
dilution of the pollutant after it has been emitted to the atmosphere to 
such a concentration that man, animals, vegetation, and materials will 
not be harmed. l 


II Control of the Pollutant at the Source: 
meee OE ene 2 Ottuvant at the Source 


Control of the pollutant at the source ma i 
y be accomplished 
keeping the pollutant from coming into existence or ohieeitye Na ae as 


existence, destroy, alter, or trap it before it can get out into the 
atmosphere, 


The most positive way to abate air pollution is + 
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pollutant from coming into existence; or if this is impossible; etea 


the quantity that does form to a bare minimum. For j 7 
in operation, this may be done es industries already 
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b) Operational changes 
c) Modification or replacement of process equipment 


d) Adopting an alternate method to produce the same 
product with fewer atmospheric wastes 


e) More efficient operation of existing equipment 


For industries developing entirely new products or processes, 
it is important to think about pollution problems in the research labora- 
tory. With such an approach: 


a) It may be found that, when looked upon from an air 
pollution angle, other lines of research differing 
from those being followed, may prove more promising. 


b) Or, by developing a more expensive process, but one 
which reduces or eliminates the need for air pollu- 
tion control equipment, it may be found that the 
overall production costs will be cheaper. 


c) If the pollutant cannot be prevented from forming, 
and its uncontrolled mass rate of emission to the 
atmosphere is excessive, equipment which destroys, 
alters, or traps the pollutant will be required. 


There are several basic methods of reducing a pollutant to 
tolerable levels of concentration before it is emitted from the stack 
top and these are destroy, mask, counteract, or collect the pollutant. 


a 


III Control by Dilution in the Atmosphere: 


The best method of controlling air pollution confines the contam- 
inant at its source. However, a second method, that of natural dilution of 
the pollutant in the atmosphere before it can reach the receptor in harmful 
concentrations, is also used. The natural dilution of contaminants can be 
accomplished by use of tall stacks, community planning in which zoning the 
use of the air isthe primary object or effect meteorological control of the 
process operations. 
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AIR POLLUTION CONTROL SECTION 





The Air Pollution Control Section is one of six sub-sections of 
the Environmental Health Services Division in the Alberta Department of 
Health and the Head of the Air Pollution Control Section reports directly 
to the Director of the Division. Besides the Head, this Section is staffed 
by four engineers, two environmental engineering technicians, four industrial 
laboratory technicians, one laboratory assistant and at least two engineer- 
ing summer students. The services of the Environmental Health Laboratory 
are also utilized by the Air Pollution Control Section. 


Basically the function of the Air Pollution Control Section is to 
develop and adminster programs for the assessment and control of air pollu- 
tion in Alberta. This is done by regulating the pollution from industrial 
plants, assessing and reducimg, where necessary, the emission of pollutants 
from various industrial operations and assessing the level of pollution in 
the vicinity of industrial operations and urban areas, Research programs 
sponsored by the Department of Health and carried out on our behalf are 
continually reviewed. These projects include the meteorological studies 
being done in Calgary and Edmonton and the mathematical modelling of dis- 
persion calculations of pollutants. 


Air Pollution Control in Alberta is governed by provincial legis- 
lation and the Regulations as such apply equally throughout the province. 


LEGISLATION AND REGULATIONS IN THE PROVINCE OF ALBERTA 


Legislation for the control of air pollution in Alberta has 
developed as the need, or anticipated need, arose. The first piece of 
legislation was to provide machinery to evaluate complaints and the remedy 
for correcting the conditions where necessary. A second item referred 
specifically to blowing dust from piles of coal and slag. These were 
enacted in 1945 and 1946 respectively. 


In 1955 the authority to make regulations was added: and in 1962 
this was amended to include pipelines as an item subject to control as to 
operation and location, 


The penalty for non-compliance with a Provincial Board of Health 
Order or with a rule or regulation for the Control of air pollution was set 
at the rate of up to $500 per day in 1955. 


The first specific regulations for the control of ai i 
r air pollution 
were formulated in 1961 and became effective on September 15th of that year 
These define air pollution and air contamination as follows: ; 


"Air Pollution" means the presence in the outdoor atmosphere 
of any air contaminant in quantities that may cause discom- 
fort to persons or endanger their health or safety, or that 


as yes 
el aia injury or damage to property or to plant or animal 
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LEGISLATION AND REGULATIONS IN THE PROVINCE OF ALBERTA: Cont'd. 


"Air Contaminant" means any solid, liquid or gas or com- 
bination of any of them in the outdoor atmosphere result- 
ing from the activities of man. 


A basic requirement of the regulations is the submission of 
plans and specifications for all new industries and incinerators (except 
those for single family residences) to the Provincial Board of Health for 
approval, Existing air pollution sources were given five years in which 
to arrange compliance with the regulations. 


The standards for air pollution control are set out in four 
sections: 
1. Smoke: as measured by Ringelmann Chart or equivalent with limits ~ 
as to density and duration in urban and rural areas. 


pores 


2. Particulates: as an emission standard of 0.85 lbs./1000 Ibs. of ¥ 
flue gas adjusted to 50% excess air for the following three ~“ 
operations: 1) Combustion, 2) Incineration, 3) Processing / 
and Manufacturing. Sd 


3. Odorous Material: as an effect standard as to use and enjoyment wee 
property and the health or safety of the public. ; 


4. Toxic or Noxious Materials: release is specifically subject to ee 
written approval by the Provincial Board of Health. Such ca 


approvals are based on the "calculated" maximum ground level a 
concentration of the material which will occur. j 


It should also be added that the Board may vary specific limits 
respecting particulates and smoke where the Board is satisfied that such 
variation is in the public interest due to location, permanency of the 
operation, specific estimated effect of the air effluent or other factors 
deemed to be pertinent. 


APPLICATION FOR APPROVAL 





The method of applying for an Air Pollution Control Approval 
varies as to detail from one industry to another, but all include a process y 


summary and flow sheet, with material balances, and information on the vs 
volume, temperature, velocity, contaminant concentration and stack height Vi 
at each point of emission. ¥ 


A follow-up to the Air Pollution Approvals is made in two ways: 


1. A periodic assessment of the emissions is made by actual measurement 
of the emissions. In the case of larger sources this is once per year. 


2, Area monitoring for the specific contaminants involved. This has been 
done in two ways up to the present, the monthly measurements such as 
exposure cylinders for H,6 and are lime papers for fluorides and 
dustfall cans for partictilates. Continuous monitoring analyzers are 
used for HoS, S05, smoke, oxides of nitrogen and oxidants. 
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ATR POLLUTION CONTROL & REQUIREMENTS AT GAS 


PROCESSING & SULFUR RECOVERY PLANTS 









Regarding air pollution at gas processing and sulfur recovery 
plants the Provincial Board of Health Regulations cover this aspect only 
generally but some o 6 acceptable pollutant concentrations resulting 
from gas plant operations are: 





1. Sulfur Dioxide: 


a)eO22 ppm as maximum ground level concentration in populated areas 
and in areas where vegetation is of commercial importance, that is, 
agricultural crops. 


b) 0.3 ppm as maximum ground level concentration in areas with little 
or relatively no population but where vegetation, such as timber, 
is potentially valuable. In areas with valuable timber the maxi- 
mum ground level concentration shall apply at tree top level in- 
stead. 


c) 0.4 ppm or higher (as conditions may warrant) as maximum ground 


level concentration in areas with no population and with scrub 
vegetation. 


2. Hydrogen Sulfide: 

Maximum ground level concentration is 0.01 ppm except in very remote 
areas. The present requirement is that H.S effluents be incinerated to 
S05 before being released to the atmosphefe, 

3. Sulfur Trioxide: 


Maximum ground level concentration is 10 micrograms per cubic meter 
(1 microgram/M’ seems achievable). 


As mentioned above the present requirement is that all tail gas from sour 
gas processing plants (gas containing hydrogen sulfide) be incinerated to 





‘convert all H5S to SO5 and then exhausted to the atmosphere at a tempera— 


ture of at least 1000°F through a stack of sufficient height so as to give 
resulting maximum ground level concentrations of SO» within the above —— 
prescribed limits, The computational methods that We use for estimating 
the expected ground level pollutant concentration and thus the estimation 
of stack height required are the Sutton Equation with the Lowry Modifica. 
tion for the case of even ground and the Pasquill Method for the case of- 
uneven ground, 5 eT ae eae : depih et 
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Another requirement is that all 
gas be flared in specially designed flare 
of fuel gas so as to ensure the successful 


Plant emergency releases of sour 
stacks with an adequate amount 
flaring of gases with low heat- 
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ing values. lxperience has shown that a minimum of 250 B.T.U./ft.> of 
gas going to the flare should be maintained. Again, depending on the 
length of flaring and the amount of gas to be flared the height of flare 

stack is usually designed so that the resulting maximum ground level 


S05 concentration does not exceed 0.2 ppm for flaring greater than one 
hour and 1 ppm for flaring less than one hour. 


Whenever condensate is stored at the gas processing plant site 
a further requirement is that the true vapor pressure of the stored con- 
densate does not exceed 12 psia. This is required so that evaporative 
hydrocarbon emissions to the atmosphere may be reduced to a minimun. / 


Another requirement is that odors be eliminated or reduced so Wi 
as not to create a nuisance in the general vicinity of the plant. R 
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APPROVAL GUIDELINES 


In order to reduce the possibility of air pollution in Alberta 
the Provincial Board of Health may require any person or corporation pro- 
posing or planning to construct a pipe line which will contain dangerous 
or noxious materials or an industrial plant, or other premises, which may 
be a source of air pollution to submit plans and specifications for the 
proposed facilities to the said Board for review before construction is 
commenced. Such plans and specifications shall be in such detail as the 
Provincial Board of Health may require, but in any event shall show the 
proposed location of the facility and its size, capacity, and the nature 
of the processing, manufacturing or other operation. No construction 
shall be undertaken or proceeded with until approval in writing of the 
proposed construction has been granted by the Provincial Board of Health 
and the Provincial Board of Health may direct such changes as to the 
location and in the plans and specifications submitted for review as the 
said Board deems necessary in the public interest. 


As regards approvals for proposed sour gas processing and sulfur 
recovery plants the procedure for approval application can be divided into 
two stages consisting of the following steps: 


Stage I - Application for Provisional Approval 


1) The company concerned notifies the Environmental Health Services 
Division of its wish to construct and operate a gas plant at a 
certain location in Alberta. The legal description and a contour 
plot plan of the proposed plant site is supplied. 


2) Preliminary plant design information, including: 


a Plant maximum raw gas inlet rate, MMSCFD (60°F & 14.7 psia). 
Raw gas hydrogen sulfide content, mole %. 

c) Plant sulfur recovery efficiency based on sulfur in raw gas, be 
d) Incinerator stack height and top internal diameter, ft. 

e) Incinerator stack flue gas volume, sulfur dioxide content, 

and temperature and velocity upon stack exit, based on plant 
maximum operating conditions. 


is prepared by the prospective owner and forwarded to the Envir- 
onmental Health Services Division. 


3) The plant owner reviews the proposed site in the field with an 
engineer from the Department of Health Air Pollution Section. 


4) The Environmental Health Services Division reviews the sub- 
mission and assesses the adequacy of the proposed facilities 
to control air pollution. Discussions may be held with com- 
pany representatives regarding questionable or contentious 


aspects. 
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Stage I - Cont'd. 


5) 


Upon the Provincial Board of Health's satisfaction that the 
preliminary design information indicates that adequate air 
pollution control measures will be instituted the Provincial 
Board of Health issues a Provisional Air Pollution Approval 
enabling the owner to proceed with the final plant design. 


Stage II - Application for Final Approval 


1) 


2) 


3) 


Final plant design information, including: 


a) Plant schematic flow sheets and plant material balances of 
the raw feed and all products and effluents. The disposal 
of all gaseous wastes shall be described in detail includ- 
ing volume, composition and temperature. Special mention 
should be made of mercaptans, hydrogen sulfide, sulfur 
dioxide and other odorous gases. 


b) Plans of incinerator stack indicating height, top internal 
diameter, location, number and size of stack sampling ports, 
sampling platform and access ladder, 110 volt electrical 
outlet location, and stack top temperature censor. location. 


c) Emergency flare stack heights and maximum amount of sour 
gases to be flared. 


is submitted to the Environmental Health Services Division. 


Other pertinent information required at this time is: 


a) A plot plan of plant area and surrounding vicinity indicating 


the Ay oe exposure cylinder (hydrogen sulfide and sulfur 
dioxide) station network, 


b) Condensate storage true vapor pressure, 


Upon the Provincial Board of Health's satisfaction that adequate 
measures will be instituted to control air pollution to acceptable 
levels for the particular plant area the Provincial Board of 
Health issues a Final Air Pollution Approval and construction and 
subsequent operation of the proposed gas processing and sulfur 
recovery plant may commence, 


The Provisional Air Pollution Approval means that the basic 


points of concern such as plant location, pollutant effluent volume, and 
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expected maximum ground level pollutant concentration have been checked 
while the Final Air Pollution Approval means that a final check on these 
basic points has been made together with a review of all the other per- 
tinent plant details and includes a list of various operating, reporting, 
and monitoring requirements once the plant is in operation. These other 
requirements cover emergency flaring of sour gases, monitoring of the 
atmosphere and stack sampling and submission of reports. 


The submission applications for approval should be sent to the 
Environmental Health Services Division, Department of Health, Administra- 
tion Building, Edmonton 6, Alberta. 


CONSIDERATIONS IN DESIGNING A SOUR GAS PROCESSING PLANT 


1) The major contributor to air pollution from a sour gas processing 
and sulfur recovery plant is the tail gas from the sulfur recovery 
plant. This gas usually has a high content of sulfur gases in the 
form of HS and S05 along with other gases such as COo, H0, 0, No 
and 05. Whenever fhe control of pollutants at their source is not 
pet re 























cool gases residue (fuel) gas is added to the incineration process 
in amounts sufficient to exhaust the incinerator flue gases to the i rel 
“atmosphere at a temperature of at least 1000°F. From the incinerator | 
the combusted tail gas is PEnPOS SONG En ry Oe ae Sloe 
to give resulting maximum ground level concentra fe) within — 
prescribed limits. The computational methods used by the Environ- 
“mental Health Services Division for the estimation of expected 
ground level pollutant concentration and thus the estimation of stack 
height required are the Sutton Equation using the Lowry Modification: \) 
“for the case of even ground and the Pasquill Method for the case of ey 
“uneven ground. 
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centrations by the Sutton Equation using the Lowry Modification 
and the Pasquill Method are included in Appendix A and copies of 
Department of Health Computer Programs No, SH505 (Calculation of || 
Ground Level Pollutant Concentrations from Elevated Sources | 
employing the Sutton Equation using the Lowry Modification) and 
No. SH510 (Calculation of Ground Level Pollutant Concentrations 
from Elevated Sources employing the Pasquill Method) are included 


in Appendix B. | ak p09 


Computational outlines for computing ground level pollutant con- | 
| 
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2) 


a) 


4) 


A second large source of air pollution at sour gas processing 

and sulfur recovery plants is the flaring of emergency releases 

of sour gas. Although this may not be considered a continual 

large source it may become one if plant upsets force the flaring 
of raw sour gas or acid gas. Since both gases have a high Hs 
content the requirement is that they be flared in specially 
designed flare stacks with an adequate amount of residue (fuel) 

gas so as to ensure the successful flaring of gases with low 
heating values. Experience has shown that a minimum of 250 B.T.U./ 
ft.? of gas going to the flare should be maintained and the flare 
pilot light should be relatively non-extinguishable so as to ensure 
that the combustible gases reaching the flare are immediately 
ignited. Depending on the length of flaring and the amount of gas 
to be flared the height of the flare stack is usually designed so 
that the resulting maximum ground level S05 concentration does not 
exceed 0.2 ppm for flaring greater than one hour and 1 ppm for 
flaring less than one hour. 


A computational example illustrating the method of flare stack 
calculation is included in Appendix C and copies of Department 

of Health Computer Rrograms No. SH508 (Calculation of Ground 

Level Pollutant Concentrations from Flares Employing the Sutton 
Equation using the Lowry Modification for the case of even ground) 
and No. SH507 (Calculation of Ground Level Pollutant Concentrations 
from Flares Employing the Pasquill Method for the case of uneven 
ground) are included in Appendix D. A graphical method for cal- 
culating the required residue (fuel) gas rates for flaring various 
rates of HS with various flare stack heights is included in 
Appendix E, 


The Sutton equation gives the maximum one-half hour ground level 
pollutant concentration to be expected on the axis of the plume 

at a distance from the source where the maximum is supposed to 
occur. Generally, the distance from the source is from ten to 
twenty effective stack heights but in any case will vary with wind 
speed and turbulence. The Pasquill Method gives the ground level 
axial pollutant concentration at any specified distance from the 
pollutant source and compensates for the difference in elevation 
between the receptor and source, if any. It is felt that the 
ground level pollutant concentration given by the Pasquill Method 
is also the maximum one-half hour concentration since good correla- 
tion is obtained when comparing with results from the Sutton 
equation but this assumption is not corroburated by the literature, 


The effective stack height is the sum of three terms - the 
actual stack height, the rise caused by the velocity of the 
stack gases, and the rise attributable to the density difference 
between the stack gases and the atmosphere, 
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CONSIDERATIONS IN DESIGNING: Cont'd. 


5) 


6) 


7) 


8) 


9) 


10) 


11) 


12) 


The effect of temperature in itself, that is, comparing winter to 
summer conditions, is not as important as the temperature gradients. 
Soon after the gases leave the stack they approach the temperature 
of the surrounding air through mixing. In the theoretical equations 
the ambient temperature is included. It has been the practice in 
Alberta to consider a temperature of 70°F. This assumption would 
make the ground concentration slightly higher than that computed 
using a winter temperature, however the summer period is the one 
which is usually considered the critical one insofar as crop damage 


' is concerned. 


The location of the maximum ground level concentration is a function 
of atmospheric stability. During unstable atmospheric conditions 
(high turbulence) the maximum occurs relatively close to the stack. 
As the atmosphere becomes increasingly stable, the location of the 
maximum ground level concentration moves further away from the stack. 


Average concentrations of a contaminant downwind from a stack are 
directly proportional to the discharge rate. An increase in dis- 
charge rate by a factor increases ground level concentrations at 
all pointsby the same factor. The discharge rate, Q, should be 
calculated at atmospheric temperature (70°F). 


Average downwind concentrations are inversely proportional to the 
square of the effective stack height. Doubling the effective 
stack height reduces maximum ground level concentrations to one 
quarter of the original value. 


For equivalent effective stack heights, average downwind con- 
centrations are inversely proportional to the wind speed. An 
increase in wind speed by a factor of two, reduces ground level 
concentrations by one-half. 


Whenever two or more stacks are closely grouped they have mutual 
influences (the plumes tend to merge into one) and the total 
pollutant ground level concentration at any point is equal to the 
sum of the pollutant ground level concentrations from the individual 
stacks at that point. 


The incinerator stack should preferably be located on the highest 
point of the plant site so as to gain as much height for pollutant 


emission as possible. 


Wherever valuable timber stands exist the allowable ground level 
pollutant concentration for that area shall apply at tree top 
level instead. 


BAG re 
CONSIDERATIONS IN DESIGNING: Cont'd. 


13) Incinerator stack sampling ports should be located about eight to 
ten diameters from any stack features which may disturb the uni- 
formity of the velocity profile, that is, bends, enlargements, 
contractions and entry or exit ports. The usual location is at about 
mid stack height. 


14) The incinerator stack sampling ports must be large enough to 
accommodate sampling probes (pitot tube, etc.) without damaging 
them when they are being inserted. The size of the sampling port 
depends largely on the thickness of the wall of the stack. For 
example a 'thin' wall may require a port of only about three inches 
diameter whereas a wall of about twelve inches may require a port 
eight inches in diameter. 


15) For small incinerator stacks (10 foot external diameter or less at 
sampling port location) usually two sampling ports at 90° to each 
other are sufficient. For stacks larger than 10 foot external 
diameter four sampling ports at 90° to each other should be provided. 
Suitable flanged port covers should also be provided. 
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COMPUTATION BY SUTTON EQUATION USING THE LOWRY MODIFICATION 


The maximum + hour downwind ground level concentration by the Sutton equa- 


tion is: ' 6 
C ppm) = 2Q 10 x 
max ol He 


= 2.35x10° (Q Ga 
tee ee 


olf 


where Cyox = maximum > hour ground concentration expected on the axis of 
the plume at a distance from the source where the maximum 
is to occur (generally ten to twenty effective stack heights 
from source), ppm by volume, 


total contaminant emission rate at 70°F, c.f.s. 

wind speed, f.p.s. 

effective stack height, ft. 

factor representing the frequency of the most frequent wind 
direction. 


FS ao 
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If the stack gases are at an elevated temperature or have an 
appreciable exit velocity from the stack, the effective stack height is 
greater than the physical stack height because of buoyancy and jetting 
action and is given by the following formula: 


He = hs +0.75 (hv + ht) 


where He = effective stack height, ft. 
hs = physical stack height, ft. 
hv = flue gas velocity rise, ft. 
ht = flue gas thermal rise, ft. 


The temperature and velocity rise are not constant for given stack con- 
ditions but depend on the wind speed and atmospheric - temperature lapse 
rate. The equation developed by Bosanquet for the maximum rise caused 
by the velocity of the stack gases is: 


ee aa ed [QT *Vs 
1+0.43 U 





U 
Vs 
where hv = maximum velocity rise, ft. 
U = wind speed, f.p.s. 
V. = exit velocity of stack flue gases, f.p.s. 
md = total stack flue gas flow rate at ambient conditions 
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SUTTON EQUATION: Cont'd. 


The Bosanquet expression for maximum thermal rise is: 


Ns] 


6,37 g QT, ¢ A e 4 
Very 


“ 


tI 


0.387 QT... hoe 
UP 


where hy; = maximum thermal rise, ft. 
T; = temperature at which stack gas density equals that of the 
ambient air (70°F), °R. 


Nem temperature difference between stack flue gases at exit 
and ambient air, “°F, 
Zeman eat 2.8 
J 
Ais ce 


aa i ac Sy da 
: oy eO.2c.¥e fae 
/QT).V; BG a 


Bre 


2 
or Mh ct t Ge It 
/ QT) -V, gG A 


G = the gradient of potential atmospheric temperature, 
EU fte 
g = acceleration due to gravity, 32.17 ft./sec.~ 


Normally only 75% of the velocity and thermal rise is used because 
utilization of the full thermal and velocity rise can never be achieved 
The critical wind speed for a given stack is estimated by calculating 
the maximum ground-level concentrations for a number of different wind 
speeds with appropriate effective stack heights. 
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The following values of constants are listed for use with the Sutton Equation 







Wind Speed m.p.h. 


u, ft./sec. 






00114 0.0041 


23-6740 





55.169 


Useful Forms 

ar = cefes. @ 70°F & 14.7 psia 

Tl = 70°F, or 530°R Vs = ft./sec. 
on, Vz = = 


1. u £t./sec. 
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COMPUTATION OF THERMAL RISE 
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COMPUTATION OF VELOCITY RISE 
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COMPUTATION BY PASQUILL METHOD 


The formula for estimating the distribution of ground con~ 
centrations downwind of a source of air~borne material as outlined in 
a paper by Pasquill is: 


Co(ppm) = 1,68 x 10° x Q 
Udoh 


contaminant emission rate at ambient conditions (70°F), cfs. 
wind velocity, f.p.S. 

distance from source along axis of plume of point of concern, 
Ct. 

= lateral spread of plume at distance 'd', degrees. 

= vertical spread of plume at distance 'd', ft. 


where 


a ao 
Hou a 


SO 
j 


This method gives the concentrations specifically in terms of the vertical 
spread of the plume (h). The magnitude of h initially increases with the 
distance d from the source and varies according to the rapidity with which 
the vertical mixing occurs in the atmosphere. Tentative estimates of h 
over a range of distance 330 ft. to 330,000 ft. are given in Figure 1, for 
6 categories of stability which are broadly specified below in terms of 
wind speed, insolation and state of sky. ‘Strong! insolation is intended 
to refer to sunny middle-of-day conditions in mid-summer, while 'slight! 
insolation refers to corresponding conditions in mid-winter. 


KEY TO STABILITY CATEGORIES IN FIGURE 1 








TABLE 1 
a Thinly 
Surface Strong Moderate | Slight Heavily Overcast <3/8 
Wind Speed, Insulation Insulation Insulation Overcast, or>4/8 cloud 
Pepe: Day or Night low cloud, night 
free night 
< 6.5 A A~B B D : = 
6.5 A-B B C D E F 
13 B BG C D D i 
20 G C-D D D D D 
> 20 GC D D D D D 
bE Rh A 5 RETIRE. | BMG 
unstable conditions neutral stable conditions 


conditions 
(For A-B take average of figures for A and B, etc.) 
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PASQUILL METHOD: Cont'd, 


The values of lateral spread @ for the two extreme distance 
s 
330 an rei ook ft. are estimated in Table 2, Values at approximately 
equal intervals between these should be taken for intermediate distan 
@,200 ft, and 33,000 ft. oi 





TABLE 2 

Stability d = 330 feet d = 330,000 feet 
Category @ Q 

A pil an 

B 40° 20° 

C 30 oe 

D 20° rOf 

E a. oh 

F 10% oo 


The axial concentration Co for an emission rate of Q c.f.s. may 
now be calculated from the Pasquill equation by substituting the distance 
d, the estimated values of h and © and the appropriate value of wind 
velocity U. The surface wind speed should be used for the short distance 
(330 ft.) and an estimated mean value between surface and geostrophic 
wind velocities throughout the main depth of the cloud for the extreme 
distance of 330,000 ft. 


If the source is elevated the concentrations at ground level 
will be reduced at the shorter distances but will tend more closely to 
the ground values as distance is increased. Correction factors Fy as 
a function of h/H, where H is the height of the source, are given in 
Table 3 below, and these should be used to miltiply the previously derived 
values of Co. For this purpose it will usually be necessary to take values 
of h at shorter intervals and this of course can be done using Figure l. 


TABLE 

h/H FV 

Gly 5.60x 107) 
0.5 +, 0 10.3 
0.6 My Ap 10_3 
0,7 0.0.2 20 
0.8 0.027 

0.9 0.057 

26 0.10 
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PASQUILL METHOD: Cont'd. 


TABLE 3 - Cont'd. 


h/H Fy 
1.5 0.36 
2.0 0.56 
3.0 0.77 
4.0 0.87 
5.0 0.91 
10.0 0.98 


EXAMPLE OF CALCULATION OF DISTRIBUTION OF GROUND CONCENTRATIONS 





Data 


Emission rate, Q= 1.15 c.f.s. of S09 
Wind velocity, U = 14.67 f.p.s. 
Stability category, C (Table 1) 
Effective stack height, H = 202 feet. 


Calculation of Co for Ground Release 


Distance in ft., d 300 3000 30,000 
Vertical Spread, h (Figure 1) L8 375 3000 
Laterial Spread, © (Table 2) 30° 30°° 30° 
Co, ppm 2065 0.39 0.0049 


Allowance for Elevation of Source 


h/H 0.5 0.8 15 2nG 4.0 
F) (Table 3) Vx Ose 000275 Ore eiuesGe 10.87 
Weert: 101 162 303 LO4 808 
ae ft. (Figs ie 690 1130 2350 3250 7100 
@ (Table ay 30° 30° 30° 30° 30° 
Co, ppm 6.3 2eh 0.62 0.33 0.077 
F) x Co, ppm 0.00063 0.065 0,223 "Os19 0.067 


Concentrations for other wind velocities are calculated in 
the same manner. 
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PROGRAM: SH505 
PROGRAM TITLE: DEPARTMENT OF HEALTH 
Stack Calculations 
Calculation of Ground Level Concentrations 
from Elevated Sources (Stacks). 
Case No. 1 — Strong Inversion 
Case No, 2 - Adiabatic Atmosphere 
DATE: December 1964 
BY: 


E. Kupchanko 


SSS ss espe he essere pt unnensnetunancsunerecanastege 


1.0 SCOPE OF PROGRAM 
1.01 Purpose 


1.011 Problem Description 


This program was designed to calculate the maximum pollu- 
tant gas concentration and the distance to the maximum 
for varying stack heights, wind velocities at any set of 
constant exit gas velocities, total gas flow rates and 
contaminant gas flow rates. It is in two parts, Case No, 
1, Moderately stable condition, and Case No. 2, moder- 
ately unstable condition. Both parts are calculated 
simultaneously. 


1.012 The advantages of this program are that it can be set up 
very quickly (2 cards for each set of conditions) and that 
it solves problems quickly. 

1,02 Method 

1.021 Method of Solution Programmed 
This air pollution program calculates the maximum pollu- 
tion gas concentration at ground level, and the distance 
from the stack at which it occurs using the Sutton equa- 
tion with the Lowry modification (1). 

Comax (ppm) = 2035x110? { Q . Am (1) 
me Hee OU 
where He = equivalent stack height defined by (2) 
He = h, + 75 (hr + ht) (2) 
where h, = stack height ifts) 
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1.0 SCOPE OF PROGRAM 


1.021 Method of Solution Programmed ~ Cont'd. 


hy = velocity rise; distance the gas rises in the 
atmosphere because of its momentum, 


hy; = thermal rise; distance the gas rises in the 
atmosphere because of its elevated temperature. 


Q = contaminant gas flow rate at ambient conditions, 
U = wind velocity (ft./sec.). 


Am = a function of wind velocity. 


Calculation of Velocity Rise, hy 
fst acchtnenaaaanbeini Rise mmLigwsit ORs SI f 


= 1 fous 
by = = ae aR dtlalls (3) 


where Qtl = total gas flow rate from stack ft.7/sec. 
at ambient conditions. 


Vs = velocity of gas in stack, ft./sec. 


Calculation of Thermal Rise, hy (Case No. 1) 
ee ae ee 


hy = 2387 Q7, A 2/3 C4) 
where A = difference T*F between air and stack gas 
temperatures. 


Z = Ind* + 2 ~-2 
sk 





— 1A 
gpl 43 /M-.28 Wm) ay 
J tl Vs eG A 


§ = acceleration due to gravity 32.17 ft./sec.~ 


G= Gm, a function of wind velocity. 


Tl = temperature at which stack gas density 
equals that of atmosphere, 
I2 = stack gas temperature, 
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1.0 SCOPE OF PROGRAM 


1.021 Method of Solution Programmed - Cont'd, 


Calculation of Thermal Rise, h, (Case No. 2) 
Serena ey Se 


Ss 
oe 
I 


= 387 Qt1A2Z/U? 


xU 
3.57VQt1-Vs 


where A = difference T*F between air and stack gas 
temperatures. 


The relation between X and Z is shown on the graph 
Bosanquet Thermal rise parameters for an adiabatic 
atmosphere, 


Calculation of Distance to Maximum Concentration - 
Cases 1 and 2 





Case 1 X 
Case 2 Xx 


20 He 
10 He (5) 


ol 


where X = distance from stack at which maximum pollutant 
concentration occurs (ft.). 


He = effective stack height. 


2.0 LIMITATIONS 


2.01 For any one set of problems the maximum number of stack heights is 
6, maximum number of wind velocities is, 10; no more than 99 sets of 
conditions may be solved at one time. Wind velocities and diffusion 
parameters are fixed although may be changed internally in the 
program. 


3.0 PROGRAM USER INFORMATION 





3.01 Input Data 








Card # Colum Format Name Unit Description 
7 pee 21 KA Number of sets of data 
1 Sad, 21 KZ Number of stack heights 
2 2-80 79H ID User's identification of 
problem 
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3,0 PROGRAM USER INFORMATION = ate 


3.01 Input Data:* Cont'd, 


Card if Column Format Name Unit Description 
B: 1-7 P7al QiL cfs Total gas flow rate at Tl 


8-14 ay ral VS fps Exit gas velocity 

15-21 F7,.0 Ie nel Temperature of gases 

22-28 F7.0 oll He Temperature at which the 
stack gas density equals 
that of atmosphere 

Pook ey | SAMA QS02 cfs Contaminant gas flow rate 


36-78 . F7.0 SE ge) fte Jth stack height (no more 
than 6) 





ID Card must always precede data card. 
3.02 Output 


3.021 Explanation of Output 





Output of Symbol Used Symbol Used in 
Input Data Units in Program Write-up 
Identification 
Stack Height ae HS(J) ie 
Exit Gas Velocity fps VS 'e 
Total Gas Flow Rate cfs QT1 281 
Wind Velocity fps U(I) U 
Stack Gas Tempera— 

ture or 2 ie 

2 

Temperature at which 
stack gas density °F Tl qT 
equals that of at- 
mosphere 
Contaminant Flow Rate cfs  QS02 Q 


3.0 PROGRAM USER INFORMATION ceKeG = 


3.02 Output 


3,021 Explanation of Output: Cont'd, 


Calculated Output Units Case 1 Case 2 Symbol Used 
in write-up 

Thermal Rise ft, ART(IL) —“BHICT) = oy 

Velocity Rise DG; AHV(I,J) AHV(I,J) H, 

Equivalent Stack ft. ES(I) ° EHS He 

Height 

Maximum PPM CMAX CMAZ Comax 

Concentration 

Distance to Maximum ft. Diets a) pis) A 

Concentration 


3.03 No programmed stops. 
4.0 OPERATION DETAILS 
4.01 No sense switch settings. 


4.02 No operator instruction card. 
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PROGRAMS: : 9H510 
PROGRAM TITLE: Department of Health 

Stack Calculations - Pasquill Method 
DATE: : August 1965 
BY: | E, Kupchanko 





1.0 SCOPE OF PROGRAM 


1.01 Purpose 
1.011 Problem Description 


This program was designed to calculate the estimation of 
downwind concentrations of air borne material. The method 
proposed gives concentrations specifically in terms of the 
vertical spread of the plume and the practical problem is 
thus transformed into making the best possible estimates of 
these quantities. 


1,012 The program uses tabular data empirically derived together 
with theoroetical equations. 


1.02 Method 


1.021 Method of Solution Programmed 
The program calculates the distribution of concentration 


downwind of a source using the method as outlined in a 
paper by Pasquill. (1) 


Co (ppm) = 1.68 x 10° x Q 
ud@h 


where Q = contaminant emission rate at ambient conditions 


u = wind velocity (fps) 

d= lateral distance from the source (ft.) 
@ = lateral spread 

h = vertical spread (ft.) 


The magnitude of h initially increases with the distance d from 
the source and varies according to the rapidity with which the 
vertical mixing occurs in the atmosphere, 


o © © K2h 
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1.022 Key to Stability Factors 


ISTAB M8 a 3 4 ) 6 
surface Strong Insolation Slight Heavily Tiiniyis< 3/8 
wind speed Moderate Overcast Overcast cloud 
fps Day or Night or>4/8 
low cloud 
Night Night 
KSTAB KSTAB KSTAB KSTAB KSTAB KSTAB 
fps 
= 6.5 1 1-2 2 Ve b a 
23 A a3 3 4 4 3) 
20 Q 3-4 4 4 4 4 
> 20 3 4 4 4 4 4 
unstable conditions neutral stable conditions 





1.023 Estimation of the Lateral Spread 9 


The values of 6 for the two extreme distances, 330 feet 
and 330,000, were extrapolated from the various observa- 
tions below. 


Stability 
Category d = 330 feet d = 330,000 feet 
(ISTAB) 
cl ee 

J. 50° pe 

- 40° 203 

3 30° ie: 

4 20% 205 

3 ig c 

6 10° D 


enn 
1.024 Allowance for Elevation of Source 


If the source is elevated the concentrations at ground 
level will be reduced at the shorter distances but will 
tend more closely to the ground source values as the 
distance is increased. A correction factor is then used 


to multiply the derived values of Co. 


a 
(2.303 H2/h@) 


e 


¢ eAaxgss 





Cont'd. 
where F = correction factor 
H = effective height of source 


= vertical spread 


1.025 Thermal Rise 
For strong, moderate, and slight insolation the thermal 
rise is obtained by the following equations for unstable 
conditions. 


hy = 0387 QtlAZ/u? 


ie 


xu 
3.57/Qt1 Vs 


whereA= difference (T°F) between air and stack gas temperatures 


QT1 = total gas flow rate from stack (ft.?/sec.) at ambient 
air conditions (70°F) 


Vs = velocity of gas in stack (ft./sec.) 


ra 
i 


distance to point of maximum pollution 


The relationship between X and 4 is shown on the graph Bosanquet 
thermal rise parameters for an adiabatic atmosphere. 


For heavy overcast day or night, thinly overcast ory4/8 low cloud, 
and<3/8 cloud night thermal rise is obtained by the following 
equations for stable conditions. 


6387 QTIA Z/u2 


I 


he 


i 


where A\ = difference (T°F) between air and stack gas temperatures 


QT1 = total gas flow rate from stack (ft.3/sec, at ambient 
air conditions (70°F) 


Z=1nJ-+2-2 
z 


V Ei: + 
ste aremesare iol [2 - you “ 4 
/@tleVs g 
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1.025 Thermal Rise 





where g = acceleration due to gravity 32.17 ft./sec. 


G= Gm, a function of wind velocity 


Tl = temperature at which stack gas density equals 
that of atmosphere (°F) 
u = wind velocity (fps) 


1.026 Velocity Rise 


We AV atone ~t tein Otl - Vs 


dot oh G/N a 
where Qtl = total flow rate from stack (ft,7/sec.) at ambient 
conditions 
Vs = velocity of gas in stack (ft./sec.) 
u = wind velocity (ft./sec.) 





2.0 LIMITATIONS 
2.01 For any one set of problems the maximum number of wind velocities 
is 5; the maximum distances from the source is ten. The two extreme 
distances from the source are 330 feet and 330,000 feet. 


3.0 PROGRAM USER INFORMATION 





3.01 Input Data 
Note: Control Card Each Set of Data, 2 - 3. 


Card # Column = Format Name Unit Description 


a 1-3 13 KA Number of sets of data 
2 1-3 ite] KWIND Number of wind speeds 
(maximum of 5) 
-6 u KDIST Number of distances 

‘ i : (maximum of 10) 
2 7-9 13 KB Lower stability number 

Se a KC Higher stability number 
2 10-12 = (KC >KB) 
2 13-18 F6.0 SH fon Height of source 


oe eo hoy 


eke] = 


3.0 PROGRAM USER INFORMATION 


3.01 Input Data: Cont'd. 

Card # Colum Format Name Unit 
‘ 2-80 79H ID 

vA 1-8 Foul QT1 efs 
4 9-16 F8.1 V fps 
4 17-24 F8,0 T2 7 
4 25-32 F8.0 ial ve 
A 33-40. F8.2 Q efs 
4 40-80 5F8.2 UL) fp 
5 1-80 10F8.0 De) oo ft. 
6 1-80 10F8,0 DECI) ft. 


Description 

Users identification of problem 
Total gas flow rate at Tl 

Exit gas velocity 

Exit temperature of the gas 
Temperature at which the gas 
density equals that of the atmo- 
sphere 

Contaminant gas flow rate 

Wind velocity (maximum of 5) 
Downwind distance (maximum of 10) 
Difference in elevations from 


corresponding downwind distance 
and base of source (maximum of 10) 





3202 Output 
3.021 Explanation of Output 
Calculated Output 


Stability. 

1-Extremely unstable 
2-Moderately unstable 
3-Slightly unstable 
4-Neutral 

5-Slightly Stable 
6-Moderately Stable 
Thermal Rise 

Velocity Rise 
Effective stack height 


Effective stack height 
at distance D(S) 


Vertical spread 
Lateral spread 


Downwind concentration 


Units 


tt. 


Symbol Used in Program 


KSTAB 
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z » 
re Pind os ia e. 

No programmed stops 
iq ! DETA LS | 









No sense switch settings 
Po: | 
No operator instruction card 


‘Program contains subroutine XXXX 
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COMPUTATIONAL EXAMPLE ILLUSTRATING METHOD OF FLARE STACK CALCULATIONS 
een se NIA VALYULALLUND 
Flaring is normally utilized during plant upset or other ab- 
normal plant condition and is a convenient means of effectively reduc- 
ing contaminant pollution to an acceptable level. The method of com- 
puting downwind pollutant concentrations from flares described in 


this paper was proposed by Professor Bouthillier and involves making 
a material and momentum balance around the flare, 


The total amount of gas going to the flare and its constituents 
are estimated, The heating value of the flared gas is calculated and if 
less than 250 B.T.U./ft.3 then additional fuel (residue) gas must be 
added to bring the heating value up to this figure. Experience has 
shown that a minimum of 250 B.T.U./ft.3 should be maintained to reduce 
the possibility of flame failure and ensure successful flaring. A com- 
plete combustion balance on all materials is then performed and the heat 
capacities of the combustion products calculated. Twenty-five percent of 
the gross available heat is assumed to be lost due to radiation. Then the 
amount of heat available for heating the excess air is calculated. The 
heat available for this purpose is the difference between the 75 percent 
of the gross heat and the heat content of the combustion products which 
at 1800°F is 35.85 B.T.U./ft.” referenced at 60°F. The total volume of 
air heated and waste gases can then be estimated. 


EXAMPLE 
The following is a detailed example of the method described 
above for the estimation of downwind ground level pollutant concentra- 


tions from flares. 


Stack Height = 65 feet. 
Stack diameter, top 6 inches, stack area 0.197 sq. ft. 


Raw Gas Entering Stack 
1,300 scfh = 384 1b./hr. 


H,8 1.20°crs 


CO,, L.207GLe 


Residue gas CH, 1.20 cfs = 4,300 scfh = 181 Tb. /hr. 


4,300 sefh = 496 1b./hr. 


12,900 scfh = 1,061 1b./hr. 


a EE: 
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EXAMPLE: Cont'd. 


Combustion Products 





2H© + 30, tench 50 oe 2505 


CO. + 0 
CH, oh 20,, D eH, 


Volumes in scfh 


4,300 HS + 6,450 Oo + 24,300 Ny, 4,300 H,0 + 4,300 SO, 
+ 24,300 N,, 

4,300 CO, 4,300 CO, 

4,300 CH, + 8,600 0, + 32,400 N, 4,300 CO, + 8,600 H.0 
+ 32,000 N, 


Total Volumes scfh 


4, 300 HAS. +74, 300 C0, Bp 818) CH + 15,000 0. + 56,700 N 


2 2 2 


4,300 SO, + 8,600 CO, + 12,900 H,0 + 56,700 N, 


Heat Released 


HS s+ ++ + 4,300 scfh x .647 M B.T.U./scf = 2,782 M B.T.U./hr. 


CH,» + «++ 4,300 scfh x 1.013 M B.T.U./sef =4,356 M B.T.U./br. 


Total 7,138 M B.T.U./hr. 


7,138,000 B.T.U./hr. 


12,900 scfh 


Average Heat Content of Gases 


553 B.T.U./scf. 


Heat Content of Combustion Products, assuming ave, of flame temp. of 1800°F. 





VOLUME ENTHALPY 
Gas sefh lb./hr. Beteue/ib. M B.T.U./hr. 
at 1800°F 
SO, 4,300 723 320 231 
ne 8,600 995 496 494 
o 12,900 610 901 550 
me 56,700) aA 170 167 1947 
82,500 6,498 ee Bye222 
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ee ee Dame = a 


. where He 


=" 3 61 


EXAMPLE: Cont'd, 


Assuming 25 percent of the gross available heat is lost due to radiation, 


Heat to excess air = 0.75 (7,138 M B.T.U./hr - 3,222 M B.T.U./hr) 


tt 


2,937,000 B.T.U./hr. 


Excess Air Warmed to 1800°F 
2 7 OOO Bets ds hr. aa 81 2 pra 
ores Bae 3925 sefh = 6,262 1b./hr. 


Velocity in Tip 


Velocity in tip = _12,900 scfh 
= 18.2 fps 
0.197 sqe ft. (3600 3 
Momentum of raw gas = Momentum of combustible products and excess air. 


Wee 1.061 1b./hr. (18.2) fos = 
s = 1,511 
6,498 1b./hr. + 6,262 1b./hr. a 


Combustion products plus excess air at 60°F 


iH 


82,500 + 81,925 = 164,425 scfh 


QT = Combustion products plus excess air at 70°F 
= 46.5 cfs 
A. = 1800° - 70° = 1730°F (stack minus ambient temperature) 


Once Qt] and Vs are obtained the downwind pollutant concentration 
can be calcalculated by using suitable thermal rise, velocity rise, and 
diffusion equations. It is considered that the thermal rise of the waste 
gases will not be as large as that resulting from the incineration of a 
similar amount of gas in an incinerator and its release from an incinerator 
stack. To compensate for this a factor of only 0.5 of the calculated thermal 
rise of the gases is used in the compution of the effective height of the 
flare stack. 


He = hs + 0.50 (hv + ht) 


= effective stack height, ft. 
= flare stack height, ft. 

hv = velocity rise, ft. 
= thermal rise, ft. 


A momentum equation is used to estimate the velocity rise. Since the volume 
of heated air and combustion products is in the order of twenty times that of 
the raw gases the velocity rise will be quite low, however the velocity of the 
gas should not be estimated at less than 1 fps. 
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2.081 Method of Solution Programmed 


A weteriel balance is mate around the flare. The nesting 
veoue of flared ges is celculated A comiete ger>astioa 
balance ot e1) wsterisle fs performed and the heet capacity 
of the combustion procucts is obtained, Twente-fiverper . 
eet of the grove syailable heat. is aesweed.to he Jost dus 
¢o rsdietion, ‘att the amount of heating excess air is 
Seloculeted, The ersileble heat for thie purpose is the 
@ifference >orwees caventy-five per cent of the gross Haat 
ani the neat comteot of the oraductse of combustion, “Tae 
Lota] volums of sir ant waste gatee can thus be obtaiowt, 
The maxitaua oolistant ge4 co 1 al ground 
ard the distsace from the stack at it occurs ip cal- 
eulated using the Sutton equation with Lowry 
(1), and the Hosanquet = Pearson method. (2), 
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PROGRAM NO. SH 508 


PROGRAM TITLE: DEPARTMENT OF HEALTH 


STACK CALCULATIONS - CALCULATION OF 


GROUND LEVEL CONCENTRATIONS FROM FLARES 


DATE: June 1965 


PROGRAMMER: E. Kupchanko 


1.0 SCOPE OF PROGRAM 


1.01 Purpose 
1.011 Problem Description 


This program was designed to calculate the maximum pollutant 
gas concentration and the distance to the maximum for varying 
stack heights, and wind velocities, using contaminate flow 
to the flare and fuel gas rates. 


1.012 The advantages of this program are that it can be set up 
quickly (3 cards for each set of conditions). 


1.02 Method 


1,021 Method of Solution Programmed 


A materiel balance is made around the flare. The heating 
velue of flared gas is calculated. A complete combustion 
balance on all materials is performed and the heat capacity 
of the combustion products is obtained. Twenty-five per 
cent of the gross available heat is assumed to be lost due 
to radiation, Next the amount of heating excess air is 
calculated, The available heat for this purpose is the 
difference between seventy-five per cent of the gross heat 
and the heat content of the products of combustion. The 
total volume of air and waste gases can thus be obtained. 
The maximum pollutant gas concentration at ground level, 
and the distance from the stack at which it occurs is cal- 
culated using the Sutton equation with the Lowry modification, 
(1), and the Bosanquet - Pearson method. (2). 


Poe ee a 
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= 5 
Co Max (ppm) = 2.35 x 10 Qa - Am (1) 
He U 
Where He = equivalent stack height defined by (2) 


He 


h, + .50 (hy + h,) 


It is not considered that the thermal rise 
obtained by flaring will be as large as that 
resulting from burning an equal amount of gas 
in an incinerator and accompanying stack. To 
compensate for this a factor of only 0.5 of 

the calculated thermal rise of the gases is 
used in the computation of the effective height 
of the flere stack, 


Where h_ = stack height (ftz) 


hy = velocity rise; distance the gas rises in the atmosphere 
because of its momentum 


h, = thermal rise; the distance the gas rises in the atmosphere 
because of its elevated temperature 


Q = contaminant gas flow rate at ambient conditions 
U = wind velocity (ft./sec. ) 
Am 


a function of wind velocity 


1.022 Calculation of Velocity Rise, hy 





4.77 ‘1 


Ppa or + hs UV, 5 


Where Qtl = total gas flow rate from stack ft.?/sec. 
at embient conditions (i.e. 70°F and 14.7 psi) 
Vs = velocity of gas in stack, ft./sec, 


© oe a eS 
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1.023 Calculation of Thermal Rise, Mt (Case No. 1) 


h, = .387 atl A 3/v? 


Where A = difference TOP, Between air and stack 
gas temperatures 


@ =inl+2-2 


J 


I= wv 3 [UL ne 28a ere. 


J/QtlVs gG emery A 





g = acceleration dye to gravity 
32.17 ft./sec. 
G = a function of wind velocity 
Tl = temperature at which stack gas equals 
that of atmosphere ( -F) 
T2 = stack gas temperature (°F) 
1.024 Calculation of Termal Rise h, (Case No, 2) 
foea ve! Qtl A 2 /U, 
Xx = xU 
3°57 7QuL ~° VS 
Where ~r = difference [°F between air and stack gas 


temperatures 


A relationship between X and Z is shown on a 
graph Bosanquet Thermal rise parameters for 
air, adiabatic atmosphere (Magill, Holden, 
Ackley, Air Pollution Handbook). 


1.025 Calculation of Distance to Maximum Concentration 
Cases 1 and 2 


Case 1 X = 20 He 
Case 2 X = 10 He 


Where X = distance from stack at which maximum 
pollutant concentration occurs 
He = effective stack height KO 


= 10. = 
2.0 LIMITATIONS 


2.01 For any one set of problems the ‘maximum number of stack 
heights is 6, maximum number of wind velocities is 10; 
no more than 99 sets of conditions may be solved at one 
time, Wind velocities and diffusion parameters are fixed 
although may be changed internally in the program. 


3.0 PROGRAM USER INFORMATION 


3.01 Input Data 


Card No, Column Format Name Unit Description 


l=2 row) KA number of sets of data 
a =k oat t KS number of stack heights 


2 - 80 79 4H ID user's indentification 
of problem 


8.0 Q quantity of gas to 
flare including fuel 
(cfh) 


3 9-14 ak iy diameter of flare (ft.) 


decimal per cent of Cl 
constituent in gas to 
flare 


6.4 C decimal per cent of C 
constituent in gas to 
flare 


3 27 = 32 6.4 C Decimal per cent of C 
constituent in gas to 
flare 


2 


3 


6.4 C decimal per cent of C 
constituent in gas to 
flare 


3 Boe Gee C decimal penecent ofec 
constituent in gas to 
flare 


4 


y) 


6.4 A decimal per cent of N 
constituent in gas to 
flare 


2 
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Card No, Column Format Name 
é pa 50 con HS 
3 eres Fe, 0, 
3 Se al ae 0, 
3 eS a co 
4, ey, ETO T, Pp 
l mee, 7.6 TF 
h tee=-56-~ 7x0. sh (9) tt. 


ID Card must always preceed data card 


Unit Description 


decimal per cent of 
H,S constituent in 
gas to flare 


Decimal per cent of 
50, contituent in 
gas to flare 


decimal per cent of 
O, constituent in 
gas to flare 


decimal per cent of 
CO, constituent in 
gas to flare 


temperature of flarg 
gases (usually 1800 F) 


temperature at which 

the density of the 
products of combustion 

at 70 F equals that of 
atmosphere (usually 70° F) 


Jth stack height 
(no more than 6) 


ae NT 2 r 


a a eb 
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PROGRAM: SH507 
PROGRAM TITLEs Department of Health 
| Flare Calculations 

Pasquill Method 
DATE: 4 April 1968 


BYs ; P, M. Ullman 


1.0. SCOPE OF PROGRAM 


1.01 Purpose 


1.011 Problem Description 


This program was designed to calculate the estimation 
of downwind concentrations of air borne material. 

The method proposed gives concentrations specifically 
in terms of the vertical spread of the plume and the 
practical problem is thus transformed into making the 
best possible estimates of these quantities. 


1.012 The program uses tabular data empirically derived 
together with theoretical equations. 
1.02 Method 
1.021 Method of Solution Programmed 
The program calculates the distribution of concentra- 
tion downwind of a source using the method as outlined 


in a paper by Pasquill. (1) 


Co (ppm) = 1.68 x 108 x Q 
‘ did Sb sh 


where Q s contaminant emission rate at ambient conditions 





u =* wind velocity (fps) 
d= lateral distance from the source (ft.) 
@* lateral spread 


h = vertical spread (ft.) 





The magnitude of h initially increases with the distance 
d from the source and varies according to the rapidity 
with which the vertical mixing occurs in the atmosphere, 


a 


eo e@ eo K 





a K47 = 


1.022 Key to Stability Factors 








ISTAB , 1 2 3 Is 5 6 
; 
Surface Heavil i 
8 y Thinly 2/5 
Wind Speed Strong Insolation Slight Overcast Overcast pa 
fps Moderate day or night or )4/8 
7 Low Cloud 
; Night Night 
i 
KSTAB KSTAB KSTAB KSTAB KSTAB KSTAB 
fps 
oo.5 1 1-2 2 rs bs 7 
6.5 1-2 2 3 4 a 6 
13 2 2-3 2 4 4 ) 
20 3 3-4 4 4 4 4 
> 20 “| 4 4 4 4 4 
unstable conditions neutral stable conditions 


1.023 Estimation of the Lateral Spread _9Q 


The values of 9 for the two extreme distances, 330 
feet and 330,000 were extrapolated from the various 
observations below. 





Stability 
Category d = 330 feet d = 330,000 feet 
(ISTAB) 

z 50° oor 

2 40° 20 

3 30; 5a 

L 20° 10° 

5 15" " 

6 LOG 5 





1.024 Allowance for Elevation of Source 





If the source is elevated the concentrations at ground 
level will be reduced at the shorter distances but will 
tend more closely to the ground source values as the 
distance is increased. A correction factor is then 
used to mltiply the derived values of Co. 


1 
F= (3,303 He/h*) 


€ 


ihe 


1.02) Allowance for Elevation of Source - Cont'd. 


where F = correction factor 
H= effective height of source 


h = vertical spread 


1.025 Thermal Rise 
For strong, moderate, and slight insolation the thermal 
rise is obtained by the following equations for unstable 
conditions. 


hy = «387 Qtl Ag / v? 


X= xu 


3257) Qt Vs 


rexswrs 


] 


where A= difference (T°F) between air and flare gas temperatures 


Qtl = total gas flow rate from flare (ft.2/sec.) at ambient 
air conditions (70°F) 


pie velocity of gas in flare (ft./sec.) 


val 
i 


distance to point of maximum pollution 

The relationship between x and % is shown on the graph Bosanquet 
thermal rise parameters for an adiabatic atmosphere. 

For heavy overcast day or night, thinly overcastor )/8 low cloud, 
and 3/8 cloud night thermal rise is obtained by the following 


equations for stable conditions. 


h= .387 QMAZ/ v2 


where A = difference (T°F) between air and flare gas temperatures 


QT1 = total gas flow rate from flare (tt.2/sec.) at ambient 
air conditions (70°F) 


Zs Ind’ +22 


J 
2 eee Eee 
J=u ( cD eG eR mere +1 
Qtl Vs 





————— 





~ KA9 


vf, 1.025 Thermal Rise - Cont'd, 


where g = acceler=tion due to gravity 32,17 ft./sec.* 


‘ G= Gm, a function of wind velocity 


Tl = temperature at which flare gas density equals 
that of atmosphere (°F) ' 


u= wind velocity (fps) 


1.026 Velocity Rise 





hv = Nat] Qtl ° Vs 
cli > o3 u/V 5 u 
where Qtl = total flow rate from flare (ft.3/sec.) at ambient 
conditions 


Vs = velocity of gas in flare (ft./sec. ) 


u = wind velocity (ft./sec.) 


1.027 Effective Stack Height 


H = hs + .50 (hv + ht) 


2.0 LIMITATIONS 


2.01 For any one set of problems the maximum number of wind velocities 
is 5; the maximum distances from the source is ten. The two 
extreme distances from the source are 330 feet and 330,000 feet, 


3.0 PROGRAM USER INFORMATION 


3.01 Input Data 
NOTE: Control Cards Each Set of Data, 2 and 3. 


Card # Column Format Name Unit Description 


af 1-3 I3 KA Number of Sets of data 
2 1-3 73 Kwind Number of wind speeds 
(maximum of 5) 
2 be I3 Kdist Number of distances 
. (maximum of 10) 
fe a 
2 sabe 3 KB Lower stability number 


ewes ks 


~ K5u 


3.0 PROGRAM USER INFORMATION 


3.01 Input Data - Cont'd. 


Card # Column Format Name Unit Description 


2 ‘10-12 13 KC Higher Stability number 
(KC KB) 
2 13-18 F6.0 SH care Height of source 
2-80 79H ID User identification of problem 
h 1-8 F8.0 Q efh Quantity and gas to flare 
includes fuel 
h 9-16 F8.4 T ite Diameter of flare 
h 17-2)  F8.0 T2 °F Exit tip temperature (1800°F) 
25-32 F8.0 T1 : Temperature at which the density of 


the products of combustion at 70°F 
equals that at atmosphere (70°F) 


h 33- 5F8.2 U(I)=I=1,5 FPS 5 wind speeds 


5 1-8 10F8.4 Cl decimal % Decimal per cent of Cl constituent 
in gas to flare 


5 9~16 C2 Decimal per cent of C2 constituent 
in gas to flare 
5 17-2) C3 Decimal per cent of C3 constituent 
| ; : in gas to flare 
5 25-32 Ch Decimal per cent of Ch constituent 
; in gas to flare 
5 33-0 cS Decimal per cent of C5 constituent 
in gas to flare 
5 41-8 Ae Decimal per cent of No constituent 
in gas to flare 
5 49-56 H,S Decimal per cent of H,S constituent 
in gas to flare 
5 57-6), SO, Decimal per cent of SO, constituent 
in gas to flare 
5 65-72 Oo Decimal per cent of O, constituent 
in gas to flare 
5 73-80 C0, Decimal per cent of CO, constituent 
y's in gas to flare & 


o.erieoe 





3.0 PROGRAM USER INFORMATION 


Card # Column Format Name 


6 
7 


h.O 


3.01 Input Data - Cont'd. 





1-80 10F8.0 DE(J) 


3.02 Output 


= ee 


Unit 
tt. 


Ite 


3.021 Explanation of Output 


Calculated Output 


Stability 
1-Extremely unstable 
2-Moderately unstable 
3-Slightly unstable 
h-Neutral 

5-Slightly stable 


6-Moderately stable 
Thermal Rise 

Velocity Rise 
Effective Stack Height 


Effective Stack Height 
at Distance D(S) 


Vertical Spread 
Lateral Spread 


Downwind Concentration 


3.03 No programmed stops 


OPERATIONAL DETATLS 


Units 


l}.01 No sense switch settings 


4.02 No operator instruction card 


Description 


Downwind distance 


Difference in Elevation 


Symbol Used in Program 


KSTAB 


AHT 
AHV 


ES(J) 


ET(J) 
VSH 
SLS 


CO(J) 


4.03 Program contains subroutine XXXX 
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FLAC JOR eSH5SO/ 
seEAG ¢ ert RAN @ SHS? 
a ee NOLON UC TOVs Mir O ys Kee sey Vi ahs file Sh SZt6 et leE SY 1) spore NE te 
WieeIRLS HIRE CECIITE apeISCS eae 
READE DS CL ACT LSS TAR) 9 ¥ CD) oS tees poles heen ee eee ie eh 


REAN( J] 9] PBIMISUSZIKSTAR SLA. SKSTAR=1 56) 9L2Z=41 97} 
READ( Ls PIETER 
12? FORMAT ngs. 


nT EE eee een nen tenon are renner — a a ee en — infeeshesgt ne 





AN US Vai 4h 
aX Cok) i 
XX(5)=14667 


xX ( L°)=44, 
ZE£01) = 4292 
WEP PE. Te 92 
Pee ee i ee she Oe ey ot a ve oe ial a eee 
PAL Gl ee 
Z£2(5)= 28930 
£2{6)=.9904 
eet? = 1 eos 
LAG) s16 7862 
£219 )=7 44964 | SE Sa ak ie a AEM Neel AO A 
27 (19) =3 61057 
DO 1990 LeleKA 
HRM, 
THOM EN, : 


ACA ey a 





anaes tchparns eo et cee np isn anos 





TU PSK. 
(© CALCULATION UF GROUND LEVEE ce 
hears az MOIN ree eee Tole FTHOD Pred G ree BAe a Bh. JLNes 1968 ote 
KWEADC]Ls 126i KWINL ski] ri SRA OK CS SHi 
i k26 FORMAT 41 F5F6 0) 
Rie ORF AD eS ee ey Rete a ROE 
| 127. FOPRMAT( AURGENs Fa.) | , 
{ Re aod SUN sd ams 84 
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~4+ 
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eof is Fee 
INCE NERA 


| S218 FITR VA 7 (FR, M gi Be fePPFRon OL Be po 

i RENE eee aa aes 9C39C49C59A2 sHI5SSO55C nae io 
102 FORMAT (1OF Req) f 
| READ CLs 119) (hh (J) 9 J1510) 

1t9 FORMAT (.10F 8.0) 

| READ CIS LIOR(DE CJ) sgets 10) 
H| baer the , 
eis. OTs a kece Gaetan SER ieee er ape sate. ol 
ViGd aes: ‘ 


| VC2=N#C? oN ieee Bi i 

















Wee ASBIKT 3 
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VOR=QHCS 
: VHZSSQ#HES i 
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are BUFO? 
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j 
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me VCOPSQRCOD | SN SAI a iri oes meee cette ar ¥ 
m VSO2SQ*SO? 
TVR sQ , 
AREHIND4e*VC1+1 792 .*VC 242 IO eRVC34IA7NHVC4EANTG Ver ie ‘3 
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i: WRITFE( 2 eAIHOC PE Ve Peed were CL — 
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I UY oD A QUAL ELNCINE RAL LON ) 
m GO) TO 996 
m CONTINUE 
me WCOV=VC1# N44 | 
PWCPHVCP# NAN? PS FERRI NINE? O° PEAR AE a” } 
i me t= 3.1194 
Wh =VC4%.15R? 
Be wCKRHEVCER 10% 
Mm WOPH=VN?X% NAA 
m WN =Vh 2s 2 DIGG 
EE a cc cee oH 
WOOZEVCO2% 61170 
mote =-VS5O2*. 1733 
oplila SEEN PEM SAW G4 WC Sell? SHUN P4+VO24WC O2eWw ane 
wee P=YC P49 #0743 eP VO 24 Ge RVC 445 HVC H4V CO? 
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INCINERATOR STACK FLUE GAS SAMPLING 


The basic purpose of conducting an incinerator stack flue gas 
sampling survey is to determine the quantity of pollutant or pollutants 
being emitted to the atmosphere at that particular time, The conditions 
of the pollutant emission are also obtained by this survey. As regards 
the Environmental Health Services Division the results of this type of 
survey are used mainly to check and see whether or not the particular 
plant is operating according to the conditions of its Provincial Board 
of Health Air Pollution Approval but other uses of such a survey are that 
it can be used by the plant operating personnel to evaluate the plant's 
sulfur recovery efficiency, to evaluate the need of air pollution control 
equipment, and to evaluate the progress of an air pollution control pro- 
gram. 


The incinerator stack sampling survey may be briefly divided 
into five separate steps which are as follows: 


1) Determining the average velocity of the flue gases at 
the sampling port location. 


2) Determining the average temperature of the flue gases 
at the sampling port location. 


3) Determining the amounts of the principal constituents 
(O05, CO, No, and C05) in the flue gases. 


4) Determining the pollutant content of the flue gases. 


5) Utilizing the above collected data to calculate the 
rate of pollutant emission to the atmosphere, 


The first step in determining the velocity of the flue gases at 
the sampling port location and thus the volume rate of flow of the flue 
ases is the selection of the velocity traverse. For circular flues 
the common case) the cross-sectional area is divided into equal areas 
laid out in concentric circles. Since the concentric circles get narrower 
as they proceed from the center to the outside of the flue, a special 
formula is used to locate the test points which are assumed to have the 
average conditions of the area they represent. The formula used is as 


follows: 


Rn = /(n-l)a + a/2 
3.1416 


the radius of the n*® point from the centre of the circle. 


kn 
aes n = n“} point from the centre of the circle. 
a = total area divided by the number of equal areas. 
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If the cross-section at the test location has a fairly uniform flow, that 
is the range of velocities does not exceed 2 to 1, the minimum number of 
traverse points is determined as follows: 


Cross-sectional Area Number of 
oguare Feet a Test Points 
Less than 2 L 
2euor25 12 
Greater than 25 20 





Where high velocity differentials or extreme turbulence of flow are en- 
countered it may be necessary to double or triple the number of points 
to establish true conditions. 


Once the flue cross-section has been properly arranged into the 
required number of equal areas the next step is to determine the average 
velocity of the flue gases at this cross-section. The formula for deter- 
mining the average velocity of the flue gases at the traverse location is: 


Vs (average) = 2.9 Fs / 29,92 x 1.00 575 x / H average 
Ps Gq 


where Vs = average velocity of the flue gas, f.p.s. 
Fs = pitot tube correction factor, 
Ps = absolute pressure of the flue gas, inches Hg. 
Gd = specific gravity of the flue gas referred to dry air at 
Ts and Ps. 
Ts = absolute temperature of the flue gas, °R. 


| fz | ave. = the average square root of the pitot tube draft gauge readings , 
taken over the entire cross-section of the flue, (inches water)2. 


The pertinent data necessary for the average velocity formula are deter- 
mined as follows: 


1. The Pitot tube velocity correction factor, (Fs). 


Two types of pitot tubes are generally used: 


A. The Standard Pitot Tube (Figure 2): 


a) This type plugs rapidly in a moist gas with a high dust con- 
centration and for this reason is limited to gases with low 
dust loadings. 


b) It has a velocity correction factor (Fs) of 0.98 to 1.02 but 


for all practical purposes it gives differential gauge read- 
ings as the true velocity head. | 
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VELOCITY DETERMINATION ASSEMBLY 
Fig. l 


STANDARD PITOT TUBE 
Fig. 2 


SPECIAL PITOT TUBE 
Fige $ 
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| factor - Cont'd. 


a) This type is used for gases which contain high dust concen- 
trations. 


b) Since it does not give a true static pressure at the point 
being measured for velocity, a correction factor of 0.83 
may be applied to theoretical velocities up to 50 f.p.s. 
without a great deal of error. However, for very accurate 
work, the tube should be calibrated over the velocity range 
of interest. 


2. Ihe Differential Draft Gauge Reading at each Traverse Point, (H). 





To indicate the velocity head measured by the pitot tube several 
types of pressure guages may be used; 


3. Ihe Absolute Temperature of the Flue Gas, 


A. 


a) Inclined draft gauge (Figures 1 and 4). 

b) U-tube. 

c) Well type manometer. 

d) Pointer type low pressure gauge (Magnehelic). 
e) Hook gauge. 


aryl 





To determine the temperature of the flue gas three instruments 
are commonly used: 


a) The standard mercury thermometer (reading above the maximum 
temperature of the gas). 

b) The standard base metal pyrometer. 

c) Ratio potentiometer. 


The temperature of the flue gas is obtained by making a temperature 
profile across the duct in the same manner as the velocity profile. 
The average of these readings is then used, 


To determine the absolute temperature, the following formula is 
used: 


Ts = 460 +T 


where Ts = absolute temperature of the flue gas, “R, 


4. Ihe Absolute Pressure of the Flue G. 





A. 


average temperature of the flue gas, ‘F. 





Ss 


To determine the pressure of the gas in the flue, two devices may 
be useds | i 


a) The flue pressure may be read on the pitot tube differential 
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4. The Absolute Pressure of the Flue Gas, (Ps): Cont'd, 





A, a) . . draft gauge by connecting the 'static' side of the pitot 


tube to the gauge, leaving the other end of the gauge open to 
the atmosphere. 


b) If the flue pressure is more than l-inch of water vacuum, a 
U-tube should be used. In using'a U-tube caré should’be taken 
to protect the open end from air currents. Such protection can 
be provided by placing a loose ball of cotton over the opening. 


B, To calculate the absolute pressure in the flue the following 
formula is used: 


Ps _almorpherse pressure + 


where Ps = absolute pressure in the flue, in. Hg. 
= gauge pressure (vacuum) in the flue, in. Hg. 


5, The Specific Gravity of the Flue Gas, (Gd). 
A. To calculate the specific gravity of the flue gas: 


a) The moisture content of the flue gas, and 
b) An analysis of its principal constituents (05, CO, No, CO.) 
on a dry basis must be determined. 


B, The moisture content can be determined by: 


a) Humidity apparatus with the proper range of thermometers. 

b) The amount of condensate collected or the volume increase 
of the bubbler solutions described in the section headed by 
'Flue Sampling Train'. 


C. The principal constituents (0,, CO, No, CO) can be determined 
by the use of an: 


a) Orsat analyzer, or 
b) Fyrite analyzers. 


Having determined the average velocity of the flue gases at 
the traverse location the volume rate of flow of the flue gases at the 


traverse location is given by: 
Q = Vs (average) x As 


where Q = volume rate of flow of flue gas, c.f.s. 
Vs (average) = average velocity of the flue gas at the traverse 


location, f.p.S. 2 
As = cross-sectional area of the flue at the traverse location, ft. 
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Incinerator Stack Flue Gas Sampling: Cont'd. 


The next step in the stack sampling survey is the determina- 
tion of the pollutant content of the flue gases. This is accomplished 
by means of a flue sampling train and for stack sampling at gas plants 
is used mainly for the collection of sulfur dioxide. 


A. The Sampling Train consists of: 


1) Sampling nozzle and sampling probe. 

2) Impingers placed in an ice bath, 

3) Silica gel colum. 

4) Dry gas meter equipped with temperature and pressure gauges on 
its outlet side. 

5) Air pump. 


A typical sampling train arrangement is shown in Figure 5. The nogzles 
and sampling probe should be made of material (stainless steel) not sub- 
ject to corrosion or damage by high temperatures. The dry gas meter 
should be protected from outside temperature changes by using insula- 
tion, or by placing it in a paper or wooden box. The rubber tubing used 
to connect the various components of the train should be heavy to avoid 
collapse due to suction within the train. It is essential that all con- 
nections be absolutely tight to prevent outside air being registered by 
the dry gas meter. A simple test for tightness is to assemble the com- 
plete train, tightly block the sampling nozzle, and carefully apply 
suction until the pressure gauge reaches its maximum. If the meter 

dial does not creep, the system is tight. 


B. Collection of Flue Gas Samples: 


1) From the sampling probe a connection is made to four gas washing 
bottles which are connected in series and placed in an ice water 
bath. The first and last gas washing bottles are fitted with 
impingers while the middle two are fitted with a coarse cylindri- 
cal frit. The first three gas washing bottles each contain approxi- 
mately 125 mls. of 5% hydrogen peroxide while the fourth is left 
dry to trap mist from the preceding three. The ice bath serves 
to reduce the temperature of the flue gas thus encouraging con- 
densation of water vapor and discouraging evaporation of bubbler 
solutions. 


2) The last gas washing bottle is connected to a silica gel column 
to reduce the humidity in the gas stream to zero for the pro- 
tection of the dry gas meter. 


3) The silica gel column is connected to the dry test meter followed 


by the vacuum pump which draws the sample through the sampling 
train, 
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FIGURE 5 


{ 


Sampling Train for Gaseous Contaminants 


|: Sampling probe 
2-lce bath 
3: Impingers 


4: Silica gel column 
5S: Thermometer 
6:Meter (dry gas) 

7 Hose clamp (Valve) 
8-Vacuum pump 
9:Pressure guage 
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B, Collection of Flue Gas Samples: Cont'd. 


4) A l-hour sample is usually withdrawn through the train at a rate 
of + cubic foot per minute for a total of 30 cubic feet. It may 
be necessary to vary the sampling rate and time depending on the 
source. 


C. Analysis of Flue Gas Samples: 


When a sample is completed the H,0, and condensed water solution is 
measured and recorded and the satiple removed to the laboratory for 
titration with NaOH. The difference between the final and original 
volumes of the sampling solution is a measure of the volume of water 
condensate. A more complete description of the laboratory procedure 
is given in Appendix A. 


After the analysis of the flue gas samples all that remains 
is to calculate amount of pollutant (sulfur dioxide) being emitted to 
the atmosphere. In addition to the flue gas velocity and volume rate 
of flow calculations already described the following calculations are 
required: 


A. Water Vapor: 


1) The volume of water vapor in the flue gas is calculated from the 
volume of condensate collected and the residual amount remaining 
in the saturated air at meter conditions (trapped in silica gel 
column) by using the following formlae: 


| V. th 
V,, = 0.00267 x Sra 
ened 
where V, = volume of water vapor (at meter conditions) condensed from 
the sample and collected in the impingers (ft./). 
V,, = volume of water collected in the impingers (cm.3). 
T, = absolute temperature of the gas at the meter (°R). 
P,} = barometric pressure (inches of Hg.). 
Mee 


m = gauge pressure (vacuum) of gas at the meter (inches of Hg.). 


2) The volume of water vapor trapped in the silica gel column at 
meter conditions is given by: 


m 
Py ate Yam 


where M, = volume of water vapor (at meter conditions ft.?). 


VP = vapor pressure of water at meter conditions (inches of Hg.). 
volume of dry gas registered by the meter (ft.-) 
barometric pressure (aes ope hea 


gauge pressure (vacuum) of gas at the meter (inches of Hg.). 


td 
son 
no tl 
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A. Water Vapor: Cont'd, 


3) The total volume of water vapor contained in the sample at meter 
conditions is then given by: 


Boe A+ eg, 
B. Sulfur Dioxide: 


1) The volume of sulfur dioxide collected by absorption in hydrogen 
peroxide converted to meter conditions is given by: 


Vigo, = AX Nx Ty x 2.4 x 107? 





ene 
where Vin, = volume of sulfur dioxide sampled (at meter conditions ft.?). 
& = volume of standard NaOH solution used (ml.). 
N = normality of NaOH solution, 
T,, = absolute temperature of the gas at the meter (°R). 
P,, = barometric pressure (inches of Hg.). 
P= gauge pressure (vacuum) of gas at the meter (inches of Hg.). 


2) The total volume of flue gases sampled at meter conditions is given 
by: 


Vi = Vn + M+ VSo5 


where Vi = total volume of stack gases sampled at meter conditions 
LT senda? «6 
volume of gas registered by the meter at meter conditions 


(Im and Pm) ° 


3) The volume percent of sulfur dioxide in the flue gas sampled is 
given by: 
VS0, x 100 


Vy 


%@ SQ» in flue gas = 


Vi, = volume of sulfur dioxide sampled (at meter conditions). 
a “2 = total volume of stack gases sampled at meter conditions 


(T,, and P,)- 
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B, Sulfur Dioxide: Cont'd. 
4) The percent moisture in the flue gas is given by: 


% moisture in flue gas = M, x 100 
V 
t 


where M, = total volume of water vapor contained in the sample 
at meter conditions. 
Vz = total volume of stack gases sampled at meter conditions 
OT and Pa) ° 


5) The weight of sulfur discharged to the atmosphere through the 
incinerator stack is given by: 

S09 

100 


Ws=Qx x 86400 x 3.7728 x 107? 


where Ws = weight of sulfur emitted to atmosphere, long tons per day. 

Q = volume rate of flow of flue gas at 60°F & 29.92 inches Hg,, 
cfs, 

SO» = percentage sulfur dioxide in the flue gas, 


INCINERATOR STACK SAMPLING REPORTS 


The results of an incinerator stack sampling survey should be 
incorporated into a report together with pertinent plant operating data 
for the time of the sampling survey. The Environmental Health Services 
Division has found it conveneient to summarize all this data in a 'plant 
operations table! an example of which is included in Table 1, The summary 
section of the stack sampling report should discuss the pollutant emission 
rate and the operation of the incinerator stack. If the plant is opera- 
ting in violation of the Provincial Board of Health Air Pollution Approval 
then the reasons for this violation and action already taken or to be 
taken to remedy the situation should be stated, 
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Ho0>_METHOD OF DETERMINING SO>_IN A SULFUR PLANT STACK 


INTRODUCTION 


This method involves bubbling a known volume of the stack gas 
through a solution of H,0,, The SO. in the gas is converted to H2S0 
which increases the Boldi cy of the sampling solution, The increase in 
acidity and hence the SO. concentration in the stack gas is measured by 


Pnenet aae an aliquot of the exposed sampling solution with standard 
aunie 


REAGENTS 


The concentration of the standard NaOH should be approximately 
O.2N. It is prepared by diluting an approximately 1N NaOH stack solution, 
The exact normality of the NaOH to be used in the titration is determined 
by titrating against standard .100N H2S0) (available in the lab) to the 
phenolphthalein end point. 


The 5% H,0, sampling solution is prepared by diluting 30% H,0o. 


Phenolphthalein indicator. 


PROCEDURE 


Before climhing the stack to make a run, approximately 125 (+ 05) ml, 
of H,0, sampling solution is measured into each of three 250 ml. impingers, 
The total volume of H,05 used mst be recorded. The impingers are placed 
in an ice-water bath just before the sampling run begins. 


After the sampling run has been completed, the impingers are taken 
into the plant laboratory where the increase in acidity of the sampling solu- 
tion is measured. There are several steps to be followed in the analysis. 


1, Measure the volume of the sampling solution, which now contains an 
unknown amount of water condensate as well as the acid, with a 
100 ml. cylinder. The difference between the final and original 
volumes of the sampling solution is a measure of the volume of water 
condensate. Transfer the sample from the impingers into a 1000 ml, 
volumetric flask. Care must be taken in rinsing the impingers and 
this rinse water added to the sample in the flask. Dilute with 
distilled water to 1 liter and mix well. 


2. Pipette an aliquot (usually 40 ml.) of the sample into a 250 ml. 
beaker containing six glass beads. 


3, Heat the solution to boiling and boil gently for ten minutes. 
Allow the sample to cool. 


Pore, ety 


~ KT 
PROCEDURE: Cont'd. 





4. Transfer the sample into a 100 ml. volumetric flask and dilute 
to the mark, Mix well. 


5, Pipette an aliquot (usually 20 ml.) into a 250 ml. beaker, add 
50 ml. of distilled water and titrate to the phenolphthalein end 
point using the standard NaOH, 


BLANK DETERMINATION 





At least one reagent blank should be run for every plant 
visited. The sample procedure that was used in analyzing the sample 
should be followed here, starting with step 2. Step 2 should be altered 
to read: 


Pipette 15 ml. of fresh H,0 into a 250 ml. beaker 
containing 6 glass beads End add 25 ml. of distilled 
water. 


The remainder of the procedure is the same. 


LAB EQUIPMENT FOR STACK SAMPLING 


One 25 ml. burette 

Two 20 ml. volumetric pipettes 

One 10 ml. volumetric pipette 

One 10 ml. graduated pipette 

One 5 ml. graduated pipette 

Two 100 ml. graduated cylinders ® 

Iwo 1000 ml. volumetric flasks with rubber stoppers 
Three 250 ml. beakers 

One 100 ml. volumetric flask with rubber stopper 
One bunsen burner with rubber connecting hose 
One retort stand 

One 4-inch cast iron support ring 

One 2~inch cast iron support ring 

One burette holder 

One 6 x 6 wire gauze square with asbestos center 
One magnetic stirrer with mouse 

One envelope of glass wool 

Two dozen glass beads 

One medicine dropper 

Stopcock grease 

Paper towels 

One 65 mm funnel 

One 500 ml. Polyethylene wash bottle 

One sheet of Kum-Kleen Labels, 
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SOURCES AND CONTROL OF POLLUTANTS ASSOCIATED WITH THE SOUR GAS AND OIL 
INDUSTRY 





I CLASSIFICATION 


A. The gaseous pollutants which are associated with the sour gas and 
oil industry may be separated into two classes: 


1. Inorganic Gases. 
2. Organic Vapors and Gases. 


B, Inorganic gas of significance as air pollutants are: 


1. Sulfur Dioxide, 
2. Hydrogen Sulfide. 


C, Organic vapors or gases of significance as air pollutants are: 


1. Hydrocarbons. 
2. Mercaptans. 


II SOURGES AND CONTROL IN PLANT AREA 





A. Major Sources 


1. Tail Gas From Sulfur Recovery 


a) The major contributor to air pollution from a sulfur plant 
is the 'tail gas from the sulfur recovery plant’. This 
gas has a relatively high content of sulfur gases ae ee 
form of HHS, $05, COS and CS,. Other gases such as ; 
HO, N, afid’0, fould be the fiain constituents of the * 


effluefit. 


b) Control Measures 





Since SO, has a higher odor threshold than H,S, the tail 
gases are burnt with fuel gas in an incinerator. In order 
to ensure relatively complete combustion of the Hos, the 
Provincial Board of Health generally requires a minimum 
stack emission temperature of 1000°F. As the tail gas 
temperature is 300 - 400°F before incineration, this ; 
higher temperature results in a significant increase in 
the thermal rise. The effluent is emitted to atmosphere 
through the sulfur plant stack, the height of which had 
been previously calculated. The height is such that at 
the approved stack operating conditions, the resulting 
calculated ground level concentration of sa? will meet 

the minimum requirements of the Provincial Board of Health. 
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II SOURCES AND GONTROL IN PLANT AREA 





A, Major Sources 


2. Hlares 


a) 


b) 


Although this is not considered to be a continual source, 
it may become one if a plant upset necessitates flaring 
of raw gas from the field or acid gas from the sulfur 
recovery unit/s. Both gases can have a high H5S content. 


Control Measures 


ace eDEED 


As in the case of the sulfur plant incinerator, the main 
purpose of the flare is to burn the HS to S0,, It is 
therefore most important that all gases reaching the flare 
tip be ignited immediately and that the flare be equipped 
with either a non-extinguishable pilot light or/and an 
automatic ignitor. In order to ensure relatively complete 
and continuous combustion of all sour streams, the Pro- 
vincial Board of Health requires that the heat content of 
any stream prior to entering the flare stack be at least 
250 B.T.U./cu. ft. It is also recommended that if sweet 
fuel is required to meet the minimum heat content require- 
ment, it be added at the base of the flare stack in order 
to ensure thorough mixing before ignition. Usually a 
liquid knock-out vessel is installed at the lower portion 
of the flare so that any liquid in the flare lines would 
be removed thus lessening the danger of extinguishing the 
flare. 


B, Minor Sources 


Although this group of sources can yield low quantities and 
concentrations of pollutants, they must be dealt with because of the 
low odor threshold of gases such as HS and mercaptans. 


1. Spent Caustic 


a) 


b) 


Caustic soda may be used to remove trace quantities of 
H,S and mercaptans from light petroleum gas, ‘The caustic 
is regenerated but some waste known as spent caustic 
remains and must be disposed of. 


Control Measures 





The spent caustic is usually incinerated in a caustic 
incinerator, Although the main product of combustion 

is sodium carbonate, there is also SO, from the trace 
sulfur gases, The combustion products are vented either 
to the main sulfur plant stack or a special stack used 
mainly for this purpose. 


oa ae 
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II SOURCES AND CONTROL IN PLANT AREA 





B, Minor Sources 


1. Spent Caustic 


b) Control Measures: Cont'd, 


The spent caustic may also be disposed of by injecting 
into a disposal well wheres proper controls will prevent 
the escape of pollutants. Care must be exercised with 
these disposal wells in order that water acquifers are 
not contaminated. Gare should also be taken to ensure 
that when the spent caustic is stored prior to injection, 
the resulting mercaptan and H,S emissions do not become 
a muisance. ‘ 


2. Spent Amine 


Plants employing an ailine process as a means of removing 
HS will have a certain amount of spent amine which, when 
discarded, will carry over HoS. The control measures are 


usually the same as for spent caustic. 
3. Burn Pits 


a) These sumps generally receive wastes which contain small 
quantities of H,S, hydrocarbons, and mercaptans dissolved 
or free. A good example of such a stream is the salt 
water from the inlet separators. 


b) Control Measures 


The best method of controlling emission from the inlet 
separator salt water or other relatively sour streams 

4s to remove the contaminates before the water is drained 
to the burn pit. As relatively large quantities of HS 
and hydrocarbons can be dissolved in the water at the 
high operating pressure of the inlet separator the water 
should be removed to a surge tank. After the major 
portion of the H.S and hydrocarbons have come out of 
solution, the date water should be stripped with either 
sweet gas or steam, 


In cases where contaminated streams cannot be treated in 
the manner outlined, a continually burning flame is 
placed over the point where the wastes flow into the burn 
pit. Any combustible gases evolved are burnt immediately 
to form S05, COos and water vapor. 


It is advisable that ground level sources of BU, ve 


minimized as funneling of the contaminant can résult in 
a high ground level concentration. 


® e ® Ls 





~ LS ~~ 


II SOURCES AND CONTROL IN PLANT AREA 


B. Minor Sources: Cont'd. 


4. Waste Pond 


a) 


Waste ponds generally receive a number of waste streams 
which for the most part carry very little in the way of 
pollutants; however, for water pollution control it is 
necessary to pond these streams rather than run them 
directly into the watercourses. Some of these waste 
streams, eg. boiler blow down water, are sent through a 
settling pond to settle out solids before entering the 
waste pond. This pond can also receive burn pit overflow. 
A number of pollutants such as H5S, S055 mercaptans and 
hydrocarbons may be in this pond. 


Control Measures 


Care should be taken to ensure that any streams which 

have a relatively high concentration of pollutants are 
treated before they flow into the waste pond as atmospheric 
dilution is the only means of controlling the pollutants 
which evolve from the pond in the form of gases or vapors. 


5. Condensate Storage 


a) 


b) 


Condensate is stored on the plant site prior to shipment. 
The main pollutant is generally hydrocarbon vapors due 

to evaporation losses of the product although if the 
condensate is not stripped, H.S or mercaptans may also 
be present. c 


Control Measures 


The Provincial Board of Health requires that either the 
stored condensates have a true vapor pressure less than or 
equal to 12 psia or if this condition cannot be met, the 
firm provide either pressurized storage or vapor recovery. 

In most cases this requirement is also dictated by economics. 
In many instances the plant can control the condensate 

vapor pressure. 


6. Pouring of Liguid Sulfur 


a) 


As the liquid sulfur leaves the sulfur recovery plant the 
major portion is poured onto a sulfur pile. The liquid 
sulfur solidifies once it has cooled, During the pouring 
and cooling sulfur vapors and hydrogen sulfide become 
airborne. The sulfur vapors result from evaporation while 
the hydrogen sulfide comes out of solution. 


os @ 19. 
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II SOURCES AND CONTROL IN PLANT AREA 


B. Minor Sources 


6. Pouring of Liquid Sulfur: Cont'd, 


b) Control Measures 


The sulfur vapor losses are dependent on the speed of the 
wind. In order to minimize losses, sulfur should be poured 
during relatively light winds. 


Although the quantities or rates of hydrogen sulfide emis- 
sions cannot be controlled, the quantities seem to be such 
that atmospheric dilution is a satisfactory control. 


7. Loading of Solid Sulfur 


a) Solid sulfur is generally loaded by mechanical shovel into 
railroad box cars for shipment. The sulfur is taken from 
the pile, dropped into a device which breaks up the large 
lumps then either conveyed or dropped into the railroad 
car. During this operation, sulfur dust becomes airborne. 


b) Control Measures 
At present, the method of control used by most firms is 
to load during only relatively light winds. In this way 


dust emissions are minimized and most of the dust will 
settle out in the vicinity of the loading operation. 


III SOURCES AND CONTROL IN FIELD OPERATIONS 





A. Well or Battery Operations 


1. Flaring 


a) Generally gas field well or batteries do not have a con- 
tinuous flare. However, well flaring may occur if the 
well requires treating or testing. 


b) Control Measures 





In the case of sour gas fields, the required flare stack 
height is based on the composition of the material to be 
flared and the expected duration and frequency of flaring. 
The general requirements are similar to those in Section I, 
A, op) except that a pilot or automatic ignitor is not 
required. Generally in such cases flaring is not permitted 


above wind speeds of 15 mph. 


cae U0 
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III SOURCES AND CONTROL IN FIELD OPERATIONS 


A. Well or Battery Operations 
2. sour Crude Storage 


a) In sour oil field areas, if the crude is relatively 
sour and it is stored on site, a vapor gathering system 
may be required. In such a system all of the tank 
vapors are gathered and sent to a flare. 


b) Control Measures 


When flaring the vapors from sour crude storage tanks 

the control measures are similar to those outlined in 
Section II, A, 2(b), except that a liquid knock out vessel 
is not required,! 


3. Salt Water Storage 


In 'wet' gas fields, salt water is removed from the gas at 
the battery or well site ina separator before the product 

is piped to the plant. As the separator generally operates 
at a high pressure, HS will be dissolved in the water. When 
the water is drained from the separator to a storage vessel 
which is at atmospheric pressure, Hs will come out of 
solution, ¥ 


b) Control Measures 





The Department of Health is presently gathering data to 
determine if the quantities of H,S evolved justify 
stringent control requirements. 


4. Burnerite 
a) In oil field areas, when lines are pigged or storage 
vessels are cleaned, the sludge is routed to the burn pit. 
for incineration, The problem is generally an aesthetic 


one resulting from the smoke emissions, 


b) Control Measures 





Section 14-3-4 of the 'Provincial Board of Health Regula- 
tions for the Control of Air Pollution! specifies the 
maximum permissible smoke emissions. However, burning 
should never take place when wind speeds are in excess of 
15 mph as the smoke can be funneled along the ground. 

If conditions for burning are ideal, the smoke should rise 
straight up and dissipate in the distance. Before burning 
a large quantity of oil field sludge, an operator should 
burn a small measurable quantity to check the dispersion, 


ve? 6 Ofe L11 


- Lil - 
IV EMERGENCY OR ACCIDENTAL RELEASES 


If any accidental or emergency releases of H5S occur, the respon- 
sible firm is to notify the Provincial Board of Health by telephone 


within 24 hours of the occurrence and in writing within 72 hours of 
the occurrence, 
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@g The air pollution survey is necessitated by an actual or potential 
air pollution problem. There must be a receptor (human, plant, animal or any 


material) which reacts or could react unfavorably to an air pollutant as well 
as a source of the pollutant. 


I TYPES OF SURVEYS 


A. Exploratory 


This is the most precise type of survey in that continuous pollutant 
readings are taken from a mobile laboratory and/or a temporary monitor- 
ing station which is located in the vicinity of the sour gas plant or 
in the sour oil field. The mobile unit or station is generally located 
at a site for a period of two to three weeks. The area to be sampled 
during any given period is generally selected on a random bases. The 
only general rule followed is that the sampling period for any one area 
be altered between summer and winter periods. The selection of the 
actual sampling site involves consideration of matters which will be 
covered in a later section. 


The only time the sampling area is not chosen randomly is if a justi- 

fiable complaint is received. Such a complaint is generally the result 

of a local citizen or group of citizens claiming vegetation damage, eye 

or nose irritation, ailing livestock or a variety of other affects of an 
aesthetic, economic, health or nuisance nature which they feel are caused 
by air pollution. 


Upon receipt of the initial complaint one or several members of the 
Air Pollution Control Section generally conduct interviews in the area of 
the complaint. If the findings indicate that the complainant had reason 
to complain and the level of pollutants in the atmosphere may be the 
cause of the complaints, then an exploratory air monitoring survey would 
be in order. 


At the present time it is felt that the following ground level con- 
centrations of S05 are the maximum which are acceptable: 
@ hour average - 0.3 ppm = 
; 1 hour average - 1.0 ppm 


3 minute average - 1.5 ppm 


The monitoring survey lasting from several weeks to over a year 
may involve one or both of the mobile laboratories as well as one or 


more temporary monitoring stations. 


B. Exposure Cylinder 


i Devartment of Health and 
This type of survey is conducted by the p 
the gas Bist or field operators in the vicinity of all sour gas plants 
or oil fields. It is a continuous type of survey in that the hydrogen 
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I Types of Surveys 
B. Exposure Cylinder: Cont'd. 


sulfide and total sulfation cylinders report what would be considered 

a monthly average pollution figure as they are exposed for approximately 

a one month period (weather permitting) before they are replaced. This 
type of survey provides the Department with a invaluable source of informa- 
tion as it is a rough indication of the trends in any given area. 


In the case of sour gas plants the Denortment attempts to set up 
the exposure stations before the plant goes on stream. In this way the 
increase in pollution over the previous ambient level can be determined. 
It is important to note that when the exposure cylinder values are studied 
in order to determine the pollution trends, only two types of valid com- 
parisons can be made. Hither each station is considered on its own or 
the yearly averages of all stations in the area are compared. As each 
area has its own unique meteorological conditions and sources and each 
eylinder is located in its own unique topography, the reported values 
are also unique. It must also be remembered that there is a general 
seasonal variance with the values being highest in the winter and lowest 
during the summer, 


There is no maximum allowable concentration as such for the total 
sulfation or hydrogen sulfide cylinders, although most of the reported 
values from the total sulfation cylinders are less than 0.3 SO, equiva- 
lent mg./day/100 sq. em. and for the hydrogen sulfide cylinders are less 
than 0.05 SC, equivalent mg./day/100 sq. em. The readings in any given 
month are cofisidered to be significant when they deviate from what is 
considered to be the norm of the station for the time considered. 

Tables I-A and I-B summarize some average readings which have been 
obtained in a few cities in Canada and Great Britain, 


As the total sulfation cylinders also react with H,5, it is advisable 
for the sake of comparing the total sulfation and hydrogen sulfide cylinder 
results, to report both results in the same units. The most common 
units are 80, equivalent mg./day/100 sq. cm. 


There is no relationship between units of SO, equivalent and parts 
per million (ppm). Although various studies have been made it is im- 
possible to relate these two unless the concentration of sulfur compounds 
is constant throughout the exposure period. This can be accomplished in 
the laboratory but it never occurs as an actual situation. As an 
example, a cylinder exposed to a low 505 concentration for one month can 
have the same final reading as a cylinder which is exposed to a relatively 
high concentration for five days. 


The exposure cylinder survey was developed in 1932 at the Building 
Research Station in England. Details of the original instrument and 
developmental investigations are described in a booklet titled 'The 
Investigation of Atmospheric Pollution!, Fuel Research Station, 1948. 
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TABLE I~A 


TOTAL SULFATION as S03 Equivalent in mg./day/dm,2 


Date Ottawa Montreal Halifax Vanacnevate 
/ Dalhousie __ Fisheries  _ B.C. _ 
f058 Jan. 0.71 2047 1.35 1.35 = 
Feb. 0.60 33 1.08 3.35 0.94 
March 0.60 2037 0.89 2.38 1.23 
April 0.52 Let? O94 2go7 z 
May 0.47 1.28 0.54 2587 0.63 
June 0,31 0.98 Oe 1.05 0.49 
daly, 0.27 0.75 O.ay 0.84 O57 
Augs (0.31 0.91 0.19 0.90 0.61 
Sept. 0.36 0.99 0.49 1.48 0.97 
Oct. 0.42 lees, 0.73 1.83 0.88 
Nov. 0.65 a0 0.49 2,66 & 
Dec. 0.83 3.18 1.03 LAT - 
1959 Jan. ras 2296 0.95 he 84 a 
Feb. .0.69 Zyl 0.67 e005 
Mar. F022 20d4 0.80 regs 
April 0.66 1.82 1,06 ores 
May 0.45 1.26 0.62 2047 
June 0.44 1.10 0.33 1.13 
July 0.30 0.80 0.10 0.67 
Aug. 0.30 0.98 0.26 0.86 
Sept. 0.29 1.04 0.20 0.94 
Oct. 0.43 237 0.48 187 
Nov. 0.68 1.69 0.94 2.09 
ine, §=860.'77 2226 1.83 2.92 
1960 Jan. 0.80 2577 0.91 PAIS 
; Feb. 0.48 2.61 1.15 239 
March 0.47 2.57 0.81 1.90 
April 0.43 1.88 0.63 1.97 
May 0.26 1.00 0.68 1.35 
June 0.30 0.71 Byes 0.81 
July 0.32 0.74 0.11 0.60 
Aug. 0.25 0.85 0.14 one 
Bape, “6.23 0.82 0.30 OFT 
Oet, 0.49 iyo O39 1.80 
Nov. 0.60 1.40 0.63 1.79 
Dec, 0459 2055 0.27 1695 
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TABLE I-B 


URBAN SITES IN GREAT BRITAIN (Data from Dept. of Scientific 


and Industrial Research 





City Year mg. S03/dm.*/day 

Summer Winter 

Cardiff City 1940-1941 0.62 Vee) 
Hall 1941-1942 Ov75 Wi? 
1942-1943 0.67 1.48 

1943-1944 0.68 TA 

Dublin 1940-1941 0.86 2.05 
1941-1942 0.58 ay 

1942-1943 0.49 0.80 

Glasgow 1941-1942 - aa 
1942-1943 i. a2 wy 

1943-1944 ~ As43 

Leeds 1939-1940 1,60 3.95 
1940-1941 1.62 3.80 

1OAI-AG LP 2.28 5.17 

1942-1943 2.67 5.62 

London 1939-1940 1.96 he 29 
(County Hall) 1940-1941 1555 5.63 
1921-TSL2 Bee Hl 5h 

1942-1943 ie Oe 4.66 

1943-1944 2.46 Le'75 

Manchester 1939-1940 1.05 42970 
1940-1941 Iie 2.87 

T9AIS1912 1.18 3.05 

1942-1943 ees 2.48 


1943-1944 0.95 3.24 
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I Types of Surveys 
B, ixposure Cylinder: Cont'd, 


This type of survey is carried out by exposing strips of lead 
dioxide and zinc acetate to the atmosphere for a known period of time, 
generally one month. Sulfur dioxide and other sulfur compounds react 
with PbO, to form PbS0, while hydrogen sulfide reacts with Zn (C,H302)9 


to form ZnS, The chemical procedures are discussed in the Sopenive 
section of this manual, 


The following facts have been established concerning the reaction 
between lead peroxide and sulfur dioxide. 


1. Both the reagent (lead peroxide) and the products (lead sulphate) 
are insoluble in water. 


2. The yield of lead sulfate is proportional to the concentration 
of sulfur diaide for all concentrations up to one part per thousand. 


3. When air containing a constant proportion of sulfur dioxide is 
blown slowly past a test surface of lead peroxide, the reaction 
goes at a constant rate until more than 15% of the lead peroxide 
has been converted. In this case the resulting reading is approxi- 
mately 13 mg. S03/100 om.</day. 


1. The reactivity is not significantly dependent on wind speed under 
the following test conditions: 


Conc,s of SO> (ppm) Wind Speed (mph) 


30-700 2067 - .67 
T=; Pee eae 
fé) js ei) mal Lies 


5, The reactivity increases 1% for a temperature increase of 12.5°F, 
6, There is no significant correlation with relative humidity. 
7, The reactivity increases when the lead peroxide surface is wet. 


8, Different batches of lead peroxide can have different reactivity. 


Although no article discussing the reaction of hydrogen sulfide and 
gine acetate could be located, the following facts were found concerning 
the reaction of lead acetate and hydrogen sulfide. As both chemicals are 
used for I,S detections it is felt that the following general statements may 


be true for each reaction: 
1. The reaction rate increases significantly with increased air turbulence. 


we oe eke 
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I Types of Surveys 


B. Exposure Cylinder: Cont'd. 


2, Variations in relative humidity between 57% and 97% have no affect on 
the rate of reaction. 


3. Methyl mercaptan, dimethyl sulfide and dimethyl disulfide in concentrations 
much higher than normally present in the outdoor atmosphere, had no signi- 
ficant effect on the rate of reaction. 


C. Research 


Sulfur_in Vegetation 


In many cases, one of the first indication that an area has an air 
pollution problem is the injury or damage to plant or animal life. 


Sulfur dioxide is the most widespread pollutant in the vicinity of 
sour gas plants or oil fields and can cause acute or chronic plant 
damage. Chronic plant injury shows up as areas of limited cell injury 
while acute plant injury results in total tissue collapse. The degree 
of plant injury caused by this pollutant follows a direct time-concen- 
tration relationship with a definite threshold. I1t appears that plants 
can tolerate a continual exposure below a threshold of about 0.2 parts 
per million without injury. It is thought that the cells in the leaf 
can in some manner detoxify such small concentrations by converting the 
sulfur dioxide into less harmful substances. However, if the threshold 

limit is exceeded, then the higher the sulfur dioxide concentrations 
the less time required for plant injury to result. 


Samples of first and second year coniferous needles and deciduous 
tree leaves are collected in plastic bags and stored in a deep freeze 
until the sulfur content is analyzed. The sulfur content as a percent 
of total dry weight of the needles or leaves is determined following 
the Official Methods of Analysis of the Association of Official Agricul- 
tural Chemists (1955). In order to include most of the growing season, 
samples are collected in the latter part of August, September or 
October. 


References: 


iL 
"Measurement of Sulfur Dioxide with the Lead Peroxide Instrument’, 
Department of Scientific and Industrial Research, 1948. 


© The Measurement of Sulphur Dioxide and Chlorides in the Atmosphere 
at Canadian Exposure Sites! - Foran, Wellington and Gibbons, CM & S, 
1958. 

3 


E. F. Galardi and Rk, M, Manganelli, Laboratory Study of Lead Acetate 
Tile Method for the Quantitative Measurement at Low Concentrations of 
Hydrogen Sulfide' - Journal of Air Pollution Control Assoc. 
133305-309 (1963). 
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II MONITORING BQUIPMENT 
A. Mobile Units 


' The Department presently has two mobile monitoring laboratories 
which are used to monitor the atmosphere in the vicinity of gas plants 
and oilfield areas, The following equipment is used in Unit No. 1: 


1. A Davis Instruments Company electroconductivity analyzer Model 
No. 7OAl for measurement of SO,. It can be operated in ranges 
of O-5 ppm and 0-50 ppm and has a lower limit of 0.005 ppm. 

A pyrolyzing furnace can be connected if requirements call for 
measuring total sulfur compounds. 


2e A Research Appliance Co, A.1I.S.1. automatic spot tape HS sampler. 


3. An anemograph MSC Type B c/w anemograph clock and drum assembly 
for indication of wind speed and direction. This equipment was 
purchased from the Department of Transport. 


4. A Science Associates Inc. Model No. 351 recording electric barograph 
for recording of barometric pressures. 


5, A Hughes Owen Co, Ltd. Model No, 8169 recorder hair hygrothermo- 
graph for recording of temperature and relative humidity. 


The equipment is the second unit is the same except for the 
following: 


1. The Davis Instruments Company unit is a Model 11-7010 RDF, 


2. The Beckman and Whitley Inc. Model K103A High Resolution Type F 
Wind system is used instead of Department of Transport equipment. 


Please note that the equipment listed above is for the sake of 
general information and does not constitute endorsement’ by this 
Department. 


III SELECTION OF SAMPLING SITE 


The choice of the actual sampling site/s will vary with the type 
- of survey to be conducted. 


For an exposure cylinder or vegetation survey, the only limitations 
on sampling sites are that they must be accessible and in the zone of pollu- 
‘tion (generally 4 to 4 miles radius from a single point source piper ie 
the quantity of SO5 emitted). As both surveys are for the purpose of indicating 
general pollution trends in all directions from the source, the entire area 
surrounding a sour gas plant or the sour oil or gas field areas themselves ; 
would be spotted with sampling sites. There may be a higher concentration o 


of eye eas 
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III SELECTION OF SAMPLING SITE: Cont'd, 


sampling sites downwind from the most predominant wind, or in a valley if it 
tended to funnel the winds. 


Sits section for a routine exploratory survey is much more in- 
volved, As tha tamporary station or mobile unit will generally be ata 
location for onl a relatively short period of time, it is most important 
that the select»i site be such that the units be in an area where the most 
information c s obtained. 

In salecting a site, some factors to be considered are: 


1, Wind Direct 


The site mus! such that the prevailing wind direction is from the source 
to the sampline site. As the prevailing winds vary with area and season 
this is a most important consideration. 

2e 4one of Pollution: 
The formulae used in the usual calculation have all been developed in 
England and ‘Inited States, While they represent a good estimation of 
actual conditions they be no means provide an absolute answer, These 
equations iniicate that the maximum ground level concentration from an 
elevated sources generally occurs from ten to twenty effective stack heights 
downwind. flowever if there is a ground level source, the maximum ground- 
level concstitration would occur closer to the source. 


3. Influence of To 





There are equations which make allowance for uneven terrain when calculating 
the resulting ground level concentrations thereby enabling an estimate of 

the resulting ground level concentration. In areas where the general terrain 
is hilly or mountainous, these equations are used, The sites where the 


resulting ground level concentration should be highest can be estimated 
from these equations. 


v Th Cx") 4s yyoe ® 
4. Exposure Cylinders: 


The results from the exposure cylinder surveys can be used to determine 
the area where the average total sulfation or hydrogen sulfide concen- 
tration is the highest. As the cylinder values are affected by all of 
the forementioned points they are a useful tool for site selection. 
However, if one intends to determine the affect of only a single source 
(eg, sulfur plant stack), care must be taken that the cylinder values 
used for site determination are not being affected by other sources 
(eg. battery flare). 


aay 2 
a3 


de Site Access 


A sampling site should be easily accessible since it must be visited 
at regular intervals. 


ee 
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III SELECTION OF SAMPLING SITE: Cont'd. 





6, Availability of Electricity: 


As most of the analytical equipment requires 110 volt power, the site 


must be serviced by electricity. Portable generators are available 
but should be used only as a last resort, 


IV PRESENTATION OF DATA 


Once an air pollution survey has begun, there is a large quantity 
of data which must be tabulated and summarized. The aim of the tabulations 
and summarization is to present the data in a form from which conclusions 
can be drawn. If the data is not presented in a form which enables the ex- 
traction of the maximum amount of information, it is not being used to its 
full capacity and is not serving the purpose of the survey. In the ensuing 
subsections, only the manner in which the Department reads and summarizes 
some of the pollution data will be presented. 


A. Exploratory Surveys 


1. Sulfur Dioxide Data 


A peak and an average reading are read from each recording on the 
strip chart. In addition to this, the date of the recording and 
the start and end time are noted. 


The method of peak reading is generally by the method of setting 
readings into triangles. Therefore, the peak reading is the value 
of the height of the triangle and the average is $ the value of 
the peak. If the reading takes the form of a rectangle the peak 


and averages are the same. 


The data is punched on computer punch cards and once the mobile 
unit has left the site or accumulated sufficient data to warrant 
analysis, all data is processed by computer. 


The computer program summarizes the average sulfur dioxide readings 
by calculating the cumulative times of occurrence in each of the 
following ranges and then calculating a frequency distribution 


table: 

S50, ppm 
0.0000 
0.0001 0.0049 
0.0050 0.0089 
0.0090 0.0099 
0.0100 0.0199 
0.0200 0.0299 
0.0300 0.0499 
0,0500 0.0799 
0.0800 0.0999 
0, 1000 0.1499 

0.5000 0.9999 eh danedaed 


1.0000 9.9950 
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IV PRESENTATION OF DATA 


A, Exploratory Surveys 


1. 


Ate 


36 


Sulfur Dioxide Data 





If the average reading is greater than .04990 ppm, information con- 
cerning the day of the reading, time and duration of reading, the 
peak and average reading, and wind information at the time of the 
reading are printed. 


All readings are correlated with winds in two distinct forms. In 
the first instance, the minutes of readings with each wind direction 
are tabulated and the percentage of the total readings which each 
wind accounted foris calculated. In the second instance, the 
number of hours of wind from each direction while the sampling 
instrument was operated is calculated and divided into the number 

of hours of readings from each wind direction. The resulting 

figure is the fraction of the winds from each direction which re- 
sulted in readings. 


In Alberta most peak SOj readings are in the 0 - 1 ppm range. There 
have been a few readings in the 1 - 10 ppm range and a very occasional 
reading in excess of 10 ppm. 


Hydrogen Sulfide Data 


The H,© data is treated in a manner similar to the SOp data. As the 
spot €ape sampler is generally run on a three hour cycle, the read- 
ings obtained are an average reading for the period. Concentrations 
observed in Alberta have ranged from zero to a level which caused 
entire tapes to blacken. 


Total Sulfation Data 


If the readings indicate that significant quantities of HS are 
present, the 505 analyzer can be equipped with a pyrolyzing furnace 
which burns all oxidizable sulfur compounds to SO,. The reported 
readings are total sulfur compounds (total sulfation). 


These readings are generally summarized in a manner similar to that 
noted for SOp and HS. 


B, Exposure Cylinder Data 


The exposure cylinders are exposed for approximately a one month period. 
When sufficient back data is available yearly and seasonal averages 

are calculated for each station and area average are calculated on a 
monthly, seasonal and yearly bases. 


The seasons are composed of the following months: 


Spring: March, April, May 

Summer: June, July, August 

Fall: September, October, November 
Winter: December, January, February 


o (ay iegueareer 
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IV PRESENTATION OF DATA 


B, Exposure Cylinder Data 


As it is difficult to detect trends from raw monthly data, the 
Department has found it advantageous to plot the cylinder results. 


At the present time plots of the monthly valuesand yearly averages 
are made for each station. 


C, Vegetation Surveys 


The data is received once per year. The values in each area are 
averaged and the yearly averages are plotted. 
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PLANT OPERATING DATA 
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PLANT OPERATING DATA 
The sulfur recovery plants and sour gas plants which have a Final 
Air Pollution Approval are required to submit to the Environmental Health 
Services Division a summary of their monthly operating data and monthly 
exposure cylinder results before the end of the month following the month 
for which the data applies. The operating information generally required 
includes the following daily averages, totals or minimums: 


1. Raw gas inlet rate (total). 

2. Percent HS in the raw gas (average). 

3. Condensate production (total), 

4. Average Reid vapor pressure of the condensate. 

5. Sulfur production (total). 

6. Sulfur or sulfur dioxide emission (total). 

7. Volume of gas flared (total for each source). 

8. Percent HoS in the flared gas (average for each source). 
9. Fuel gas added to flared gas (total for each source). 
10. Stack emission temperature (minimum). 


When the data is received, the daily sulfur emission rate is calculated 
according to the following general equation: 


Sulfur emission rate = (Raw gas inlet x % HAs x conversion factor) - 
LTS MMSCFD production 
LTS 


The conversion factor is calculated as follows: 


a ee ae yk = 3.7728 x 107°? _LTS 
1b. mole H25) 4n9 13 ¢t.3 HoS 2240 lbs. S tie HS 
sole Hos 


The ideal gas law is used to adjust the conversion factor to the applicable 
raw gas inlet conditions. 


From the submitted monthly operating data, the Division plots 
the following information on running monthly - yearly graph paper 5 Note 
that any figures involving a 'Percent of Time! are calculated by dividing 
the number of days the condition was not met by the total number of operating 


days: 


1. Average sulfur dioxide or sulfur emission rate: 
a) Submitted emission rate 
b) Rate calculated by Division 
2. Percent of time emission limit exceeded: 
a) Submitted emission rate 
b) Rate calculated by Division 
36 aie of time emission temperature was even momentarily below approval 
minimum. Prien 
4. Average of the minimum emission temperatures. 


_ 
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APPENDIX A 





LEAD DIOXIDE EXPOSURE METHOD FOR THE DETERMINATION OF TOTAL SULFATION 
nt ae AS SE NALLON OF TOTAL SULFATION 


In exposure tests for the determination of mean values of total 
sulfation pollution in the atmosphere the following steps are carried out: 


1) Strips of surgical gauze two feet long and two inches wide 
are wrapped around & oz. Armstrong jars and secured to the 
Jars with cotton thread, 


2) The lead dioxide paste is prepared by mixing approxi- 
mately 100 gm. of PbO, powder with a solution of gum 
tragacanthin., The oie of gum tragacanthin is pre- 
pared by mixing 3 gm. of gum tragacanthin with 15 mls. 
methanol and diluted to 300 mls. with distilled water. 
This gum soltuion is added to and mixed with the PbO5 
powder until a consistency is reached that affords easy 
spreading of the PbO}, paste on the surgical gauze with 
a small brush. If approximately 100 gm. PbO» are used 
to make the paste of the desired consistency, there 
should be enough paste to cover 10 strips of surgical 
£auZee 


GRAVIMETRIC ANALYSIS FOR TOTAL SULFATION FROM EXPOSURE CYLINDERS 





The gauze strip from the jar is removed and placed in a 600 ml. 
beaker after the exposed surface area has been measured. Care must be 
exercised in removing the gauze strip and transferring it to the 600 ml, 
beaker or some of the PbO, may come off in this process. Add 5 gm. of 
NagCO03z to each sample foltowed by 100 ml. of distilled water. The mixture 
is stirred and allowed to stand for at least three hours, and then heated 
at 80°C for between 30 and 45 minutes, without letting the volume to go 
below 80 ml. The mixture is then filtered with Watman No. 12, 18.5 cm. 
folded filter paper doing an adequate job. Wash the contents of the 
filter paper several times with hot distilled water and acidify with 
1:1 HCl using methyl orange as an indicator and while stirring add about 
20 ml. of warm 10% BaCl5 solution. The precipitate is digested at 
80 - 90°C for at least two hours and then filtered through Watman No, 40 
ashless filter paper. After the precipitate is washed with small portions 
of hot water to make the washings chlorine free, the filter paper and 
precipitate are placed in weighed crucibles that have been pre-ignited 
at 800°C for at least 3 minutes. The filter paper is charred by the use 
of a bunsen burner and ignited for not less than thirty minutes at 
800°C in a muffle furnace. Cool the crucibles and weigh to determine 


the weight of the precipitate. 
Peary ne) 


a ee 


CALCULATION 


2 Shh +d 80,06 x 100.0 
S03 equivalent mg. /day/100 sq. cm. = mg. of Ba pay x 253 )0 


days x area 
In the above formula 
'days' represents the number of days that the cylinder was exposed. 


‘area! represents the exposed surface area of the gauze strip. 


ye represents mg, of SOQ, in the right hand side of the equation. 


'100_' 


re represents the ‘per 100 sq. cm.' in the right hand side of the 
ea 


equation. 
'80,.06' represents equivalent weight of S03- 


1233.42' represents equivalent weight of BaSO) e 


« we)  T, 


wi? oe 
ZINC ACETATE EXPOSURE METHOD FOR THE DETERMINATION OF HYDROGEN SULFIDE 
—————— ee SE NATION OF HYDROGEN SULFIDE 


In the exposure tests for the determination of the mean value of 
hydrogen sulfide pollution in the atmosphere, two inch wide strips of Watman 
No. 2 filter paper are impregnated with a solution containing 5% zinc acetate 
and 5% glycerine. The filter paper is soaked for fifteen to twenty hours in 
the acetate solution and then removed and allowed to dry until the filter 


paper strips are just slightly damp. The tapes are then secured to the jars 
with cotton thread. 


VOLUMETRIC ANALYSIS OF HYDROGEN SULFIDE FROM EXPOSURE CYLINDERS 





The exposed surface area of the paper is measured before it is 
removed from the jar, After the paper is removed it is folded to fit into 
the bottom of a 600 ml. beaker. 150 ml, of distilled water is added thor- 
oughly immersing the filter paper, then 10 ml. of 0.025 N iodine followed 

by 5 ml. of concentrated HCl, If, after stirring the brown color of the 
solution fades out, 10 ml. portions of iodine are immediately added until 
the brown color remains and record the total amount. The solution is back 
titrated with 0.025N sodium thiosulfate adding + ml. of starch solution 
near the end point. The end point is blue to clear to pink, The titration 
volumes should be read and recorded to the nearest 1/10 ml. An iodine 
determination should be made on the blank, 


All volumes of iodine are converted to equivalent volumes of 
Na So03 eg. if 9.7 ml. of 0.025N NaoS,03 are found to equal 10 ml. of 
iodine (9.7/10 and if 20 ml. of iodine’wére used initially, this would 
be recorded as (20 x 9.7/10) 19.4 ml. of NayS2036 


Standardization of 0.025N Sodium Thiosulfate 


Measure 10 ml. of 0.1N potassium dichromate solution into a 
400 ml. beaker containing 25 ml. H,0, 5 ml. concentrated HCl and 1g gm, KI. 
Allow the reaction to proceed in the dark for 5 minutes, then dilute to 
500 ml. and titrate with 0.1N NajS,03 adding ¢ ml. starch solution near 
the end point. The end point changes from opaque blue to a transparent 
pale green. 


Normality of NajS.03 = (ml, of potassium dichromate) x (normality of 


potassium dichromate) 


CALCULATIONS : 
1. mg. of S03/100 sq. cm. = ml, x N x 40,03 x 100 


area in sq. cm. 


In the equation 


'ml.! represents the number of mls. of iodine that reacted with the 
sulfite in the sample expressed in terms of ml. of 0,025N : 
NapS.03. For example, if in running the iodine blank you fin 


oes 
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CALCULATIONS: Cont'd. 


tml.' . . that 9.7 ml. of 0.025N Na.S,0, are equivalent to 10 ml. 
of the iodine which was Aged Sor aaah sample. If after 

adding 10 ml. of iodine to the sample it is found that it 

takes 4.7 ml. of 0.025N Na 550 to back titrate the excess 

iodine so in this case the sample contained the sulfide 

equivalent to (9.7 - 4.7) or 5.0 ml. of 0,025N NagS,03 

and in this case the 'mi.' would be 5.0 


'N®? represents the normality of the NajSo03. 
'40.03' represents the equivalent weight of S03. 


‘ml, x N x 40.03' represents the 'mg. of S03! in the right hand side 
of the equation, 


2. Corrected mg. of SO, /100 sq. cm. = mg. of S0,/100 sq. cm. (of sample). 
~mg. of S03/100 sq. cm. (of blank). 


3, Mes. S03/day/100 sq.e cm = corrected mg, of 503/100 sq» cm. 
no. of days sample was exposed. 











